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I. 



INTRODUCTION 



This study was made to present a method for the conceptual 
design of tactical missiles. The starting point for the design 
was a recently completed report by General Dynamics, Convair 
Division entitled Rapid Approach to Missile Synthesis (RAMS) . A 
procedure was then developed more akin to aircraft design and 
which bears little resemblance to RAMS and which uses basic 
equations to size components rather than nomograms and table 
look-up. The procedure starts with a threat description and 
proceeds with the formulation of performance objectives and 
a conceptual design of a tactical missile to counter the threat 
The design is not unique, but as will be shown, is a compromise 
of parameters to give one possible solution to the design prob 
lem; therefore, the point design arrived at is not necessarily 
the optimum design for the presented threat. An attempt is 
made to find the optimum performance within specific areas of 
the design process. 

Throughout this study the theory involved is explained and 
specific examples are worked. A complete design example is 
worked out in detail. It is an air-to-air missile designed to 
counter the new Soviet RAM-K fighter aircraft. This example 
is worked in each section immediately following the develop- 
ment of the theory for that specific area. 
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II. PROBLEM DEFINITION 



A. THREAT ANALYSIS 

1. Operational Requirements 

The design of a new missile is usually in response to 
an operational requirement which arises as the result of one 
of the following: (1) A new technology provides the means to 

design a more effective missile to meet a current threat. An 
example of this might be an advance in material science, which 
allows higher inlet turbine temperatures for a turbojet engine; 
therefore, allowing higher missile flight speeds. (2) Intelli- 
gence indicates a new threat for which existing missiles are 
not effective. (3) Operational reports indicate a current mis- 
sile is inadequate against a current threat. 

Regardless of how the operational requirement is derived, 
a statement of the threat is required before the design process 
can proceed. Experience has proven that one missile cannot be 
designed to meet all types of threats without seriously com- 
promising performance or effectiveness. This can be illustrated 
with the design of the warhead. A contact fuze, shaped charge 
warhead designed to penetrate and kill hard targets such as 
tanks, would not be effective against a highly maneuverable 
aircraft for which the expected miss distance is several feet. 

For this reason the design of a missile must start with a de- 
tailed analysis of the threat. The more detailed this analysis 
is, the more effective the final design can be. 
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2 . Design Example (Operational Requirement) 

A design example will be used as a continuous thread 
throughout this thesis to demonstrate applications of the theory. 
An air-to-air missile will be designed to counter the new class 
of Soviet fighter, which is in the advanced development stage 
at the Ramenskoye Experimentation Center. The fighter, as des- 
cribed in Aviation Week[l,2]has been designated the RAM-K. The 
RAM-K is a twin engine fighter with variable geometry inlets and 
swing wings. The aircraft bears a resemblance to both the F-14 
and F-15. It is expected to be the recipient of a new look 
down, shoot down radar and the 40 km range AA-X-9 missile. The 
following unclassified dimensions and performance data are avail- 
able on the RAM-K: 

Wing span 
Overall length 
Gross weight 
Maximum speed 
Service ceiling 
Figure (2-1) is a drawing of the RAM-K. 

3 . Scenario 

The scenario within which the missile is expected to 
operate should also be described. If the normal mode of operation 
of the threat is not known, an attempt should be made to define 
the most demanding scenario that can be expected. For a defen- 
sive weapon the most challenging incoming threat will normally 
be a head-on encounter. The threat profile may vary from a high 
level attack with a terminal dive to a low level attack with a 



40 ft 
64 ft 
6Q000 lbs 
M=2 . 5 
60/) 00 ft 
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RAM-K 



} 



Figure 



( 2 - 1 ) 



ram-K Fighter. 
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terminal pop-up maneuver. In the case of the AS-6 (Kingfish) 
two modes of attack can be expected. In such a case both pro- 
files must be evaluated to determine the most demanding in 
terms of missile performance objectives. For an offensive sys- 
tem, such as an air-to-air missile designed to intercept and 
destroy an enemy fighter before it launches its weapons, the 
scenarios analyzed should include all possible encounter 
geometries . 

4 . Design Example (Scenario) 

The scenario for the above threat would likely be an 
intercept situation in defense of the fleet high value unit. 

The scenario is taken to be a head-on encounter with the missile 
and the target at the same altitude. Since the combat specifi- 
cations of an aircraft are normally given at 10,000 ft., this is 
taken as the scenario altitude. 

B. HISTORICAL SURVEY 

Missile design is an iterative process, and the first time 
through the design loop many assumptions have to be made con- 
cerning component sizes and weights. One method of approach 
at the early stage is to employ historical data of existing mis- 
sile sizes and weights; since justifications for these parameters 
were made duirng their design processes. An example of the use 
of historical data in determining the initial missile length 
can be made with the length to diameter ratio. The length to 
diameter ratios of existing missiles of the same class as that 
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being designed are collected, and an average is computed. The 
diameter of the design is fixed by one of three driving factors 
(propulsion, warhead, or guidance) . From the average length to 
diameter ratio the initial missile length is then estimated. 

From this historical data, initial choices based on the experi- 
ence of others can be made for many of the missile parameters. 
These parameters define a baseline missile, which is the initial 
configuration from which design iterations and refinements can 
be made. 

Since missiles are designed for specific missions, speci- 
fic parameters such as length and diameter are of little value 
in comparing missiles. Dimensionless ratios such as length to 
diameter, L/D, ratio and aspect ratio, AR, are more meaningful 
when relating missiles. Some parameters which are useful in 
defining the baseline missile are listed below: 

L/D = Length to diameter ratio 

L^/D = Nose length to diameter ratio 

AR^ = Aspect ratio of the wing 

AR^ = Aspect ratio of the tail 

W/S = Weight to lifting surface area ratio 

V. = S^l V (S ,-d -) Tail volume coefficient 

t t t ref ref 

W^W^ = Gross weight to warhead weight ratio 
The tail volume coefficient, V t , is a dimensionless para- 
meter used to initially size the tail. For a tail control mis- 
sile, it is a measure of the relative control effectiveness when 
comparing missiles. For a wing control missile it is a relative 
measure of stability. 
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A complete historical survey should not be limited to 
the parameters listed here. Any dimensionless parameter which 
will add information about the proposed design should be inclu- 
ded for completeness. 

Table 2-1 is an example of a collection of such para- 
meters for existing air-to-air missiles [1, 3] . in this table 
the subscript, c, is used to indicate a canard control surface. 

C. LAUNCH PLATFORMS AND PHYSICAL CONSTRAINTS 

The problem definition must also include a description of 
the intended launch platform for the missile. The aircraft or 
shipboard launcher from which the missile will be launched will 
fix many design features of the missile. For instance the most 
important consideration in the problem definition phase is any 
physical constraints imposed by the launcher. Since it is not 
normally feasible, economically to design a launcher to fit the 
missile, most new missiles must fit existing launchers. In the 
case of shipboard launchers there will be a maximum length and 
diameter and a maximum launch weight which can be accomodated. 

For the case of an air launched missile, there will be a 
maximum weight, and the dimensions may also be limited due to 
the performance requirements of the aircraft. 

For the air-to-air missile design example of this study, 
the launch platforms will be the F-16 and F^lSA. Figure (2-2) , 
which is from Interavia [4, 5]shows the pylon weight limitations 
for these aircraft. From these figures it can be seen that the 
maximum launch weight for this design is limited to 2500 pounds 
by the pylon limitations of the F-18A. 
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Load Carrying Capability of F-18A 




Load Carrying Capability of F-16 




3500/5.5 250/9.0 

2200/5.5 2500/5.5 



Figure (2-2). Weight Limitations. 
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D. MISSION PROFILES AND PERFORMANCE OBJECTIVES 



The mission profile of a missile consists of dividing the 
flight into fundamental segments which consist of a single func- 
tion, such as boost to cruise speed and altitude, cruise to 
the target, and terminal homing phases. The mission profile 
will vary from missile to missile. For a cruise missile it 
may consist of a series of pop-up maneuvers and low level cruises. 
For a short range missile, the profile may be entirely terminal 
homing. The Mach number and altitude are specified at the begin- 
ning and end of each mission segment, as well as the range 
covered by each segment. The range covered by a segment can be 
considered in one of two ways. If the segment is short, such 
as the terminal phase, the distance along the intended flight 
path is considered. For longer range missiles, the distance over 
the ground is of importance. The mission profile must be defined 
during the problem definition phase in order to specify missile 
performance objectives. 

The mission profile of a missile normally consists of a 
boost, cruise (mid-course) and terminal phase. The boost phase 
accelerates the missile to its flight speed. This acceleration 
may be large for a surface-to-air missile which must be acceler- 
ated from rest to a high supersonic speed or it may be small for 
an air launched missile which has the speed advantage of the air- 
craft from which it is launched. 

The cruise segment, or mid-course phase, primarily is used 
to deliver the missile to a point in space where the seeker can 
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acquire the target. The range and speed of the cruise segment 
is then a function of the required stand-off distance for the 
target. 

The terminal phase is somewhat more difficult to analyze. 

If the target maintains constant heading and velocity the flight 
path of the missile may be modeled by circular arc segments. This 
method will give approximate values of range. If the target 
maneuvers, the missile must follow and the range and speed require- 
ments become more complicated. 

In determining the range and speed requirements, all expected 
encounter geometries should be analyzed. The most demanding 
encounter will then fix the performance objectives. The most 
demanding speed requirement, in terms of maintaining a minimum 
stand-off distance, will normally be a head-on encounter. Al- 
though the required missile speed and range are determined in 
this section, the missile velocity may be varied later in the 
design process due to guidance considerations. 

1 . Design Example (Mission Profile) 

From the threat defined above, the ideal situation would 
be to obtain a fire control solution and launch such that the 
minimum separation distance between the launch aircraft and the 
target is 40 km. This can be accomplished in one of two ways. 

The launch aircraft can fire a semi-active homing missile at 
such a range and speed that intercept occurs before the minimum 
range is reached, or an active homing missile can be fired, and 
once missile lock-on is achieved, the launch aircraft can maneu- 
ver to maintain the minimum separation distance. 
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The active radar homing missile would decrease the 
launch range and the range required of the missile, which will 
lessen the cons taints on the launch aircraft. The terminal por- 
tion of the engagement is a function of the guidance law and will 
be determined in Chapter 3. 

Two cases are investigated to determine the effect on 
the range requirement of the missile when a minimum separation 
distance from the target to the missile of 40 km is maintained. 
The first case is a semi-active homing missile, for which the 
launch aircraft must maintain a closing course until intercept. 
The second case is an active homing missile which has a lock-on 
range of 10 km. The launch aircraft may then maneuver to main- 
tain a separation distance. 

a. Case 1: Serai-active homing missile 




< 



‘ 7 

Rq = Range at which the missile is launched 

M^ =1.5= Launch Mach number 

M._ = 2.5 = Missile Cruise Mach number 
M 

a = Speed of sound 

V = M_ a = 1.5a = Launch speed 
Li Li 

V T = M^a = 2.5a = Target speed 

V._ = M „, a = 2.5a = Missile speed 
MM 
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The instantaneous range from the missile to the target is 



given by, R^ 



KmT - R 0 - (V M + V 1 



( 1 ) 



The instantaneous range from launch aircraft to the target 
is given by, R LT 



"LT - R 0 - <V L + V 



( 2 ) 



If the target does not maneuver, intercept will occur at 
t f , when = 0 



0 = 



t, = 



R 0 - (V M + 



_ 



<W 5a 



If the launch aircraft is at the minimum separation distance, 
R^ t = 40 km, when intercept occurs. 



40 km = R Q - (V L + V T )t f 



substituting for V^/V^ an< ^ **f 



R 



40 km = R 0 - (1.5 + 2.5)a(^) 



Solving for R 



Rq = 200 km = Launch range 
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The range required of the missile is then, R^. 

= V^t^ = 100 km = 53.96 nmiles 

If the missile speed is increased to M__ = 3.0, then 

M 

R q = 148.15 k and, 

R^ = 80.81 km = 43.61 nmiles. 

b. Case 2: Active Homing Missile 

The lock-on range is a function of the seeker in the 
missile and will be covered later in this thesis. If it is 
assumed that the launch aircraft must maintain its course until 
lock-on occurs at a range of R^, the problem can still be solved. 
The geometry is the same as in Case 1. Instead of following a 
constant course until intercept, the launch aircraft must now 
only maintain a closing course until R^ = R^. Then from 
equation (1) 



'W - "LO = R 0 - (V M + V'fl 

Solving for t^ 

fl <VV 

If at the time of target lock-on, t^, the target 
and launch aircraft are at the minimum separation distance, 
Rlt = R m in' from equation (2) , 
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R min R 0 



(V T 



V T }t fl 



Inserting for t^, 



R min R 0 



- <VV 



(r o- r lo )/( V v t> 




(3) 



For the same geometry and relative speeds of the 

first case, with M„, = 3.0, 

M 

t = 

fl 5.5a 

A reasonable value of lock-on range is 10 km. This 
will be shown later in the guidance section of the study. The 
time to lock-on then becomes , 

R -10 

t = —2 

fl 5.5a 

From equation (3), Rq then becomes, R Q = 118.52 km . The missile 
range to lock-on is then, R^, 

R M . = V t... = 59.19 km . 

Ml M f 1 
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If the target does not maneuver the time from lock- 



on to intercept becomes, t^* 



t = !i 0 - 

« <v M+v 

and the missile range from lock-on to intercept becomes, 

Rm 2 = V M fc f2 ~ 5.45 km 

The total missile range is then the sum of the two, 
R^ = + R M 2 = 64.64 km = 34.88 nmiles . 

As can be seen from the above analysis, both the 
detection range of the target and the required missile range are 
decreased significantly when an active homing missile is used. 

On the other hand, it must also be remembered that the complexity 
and cost of the missile will be increased as a result of choosing 
an active radar seeker. For the design example in this study, 
an active radar seeker is chosen; therefore, the maximum range 
requirement will be 35 nmiles at a speed of = 3.0, however, 
this missile velocity is tentative until a guidance analysis 
is complete. 

From the preceding analysis the mission profile is 
determined. It must be kept in mind that the mission profile 
may be changed during the design process to meet other design 
objectives. The following profile assumes both the target and 
launch aircraft at the same altitude. 
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BOOST 



CRUISE 



TERMINAL 

HOMING 



A 






+ 



if 





SEGMENT 


M. 

begin 


h 

begin 


M ^ 
end 


h , 
end 




RANGE 


(1) 


Boost 


1.5 


10,000 ft 


3.0 


10,000 


ft 




(2) 


Cruise 


3.0 


10,000 ft 


3.0 


10,000 


ft 


29.6 nmiles 


(3) 


Terminal 


3.0 


10,000 ft 


3.0 


10,000 


ft 


5.4 nmiles 
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III. GUIDANCE LAW SELECTION 



Although the specifics of the guidance system is beyond the 
scope of conceptual design, the selection of a guidance law is 
necessary for initial calculations. The warhead design depends 
on the expected miss distance between the missile and the target 
and the lifting surface area depends on the maneuvering require- 
ments of the missile. Both the miss distance and maximum acceler- 
ation required are functions of the missile guidance law. 

The guidance law for a missile is the analytical formulation 
used by the guidance system to convert sensed target information 
into missile steering commands. Three general guidance laws are 
used. Most others can be forced to fit into one of these cate- 
gories. These are: 

1) Pursuit Guidance 

2) Line-of-Sight Guidance 

3) Proportional Guidance 

A. PURSUIT GUIDANCE 

A pursuit guidance law is illustrated in Figure (3-la) and 
is one in which the missile velocity vector is always directed 
toward the target. The target and the missile velocity vectors 
must therefore be sensed; so this type of guidance normally 
assumes an on-board tracker. The missile may have a separate 
mid-course guidance package to increase range, but target lock- 
on initiates the pursuit guidance for the terminal homing phase. 



35 




Figure (3-la). Pursuit guidance. 



X 




Figure (3-lb) . Beam rider guidance. 





Figure (3-lc) . Proportional guidance. 
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For this reason it has the advantage of launch-and-forget at 
lock-on. Since the signal processing is limited to looking and 
pointing, the avionics are relatively simple and usually on- 
board the missile. An option for this type of guidance would 
be to include a lead angle to accomodate faster moving targets. 

B. LINE-OF-SIGHT GUIDANCE 

Line-of-sight guidance is used in a beam rider missile. This 
guidance scheme is illustrated in Figure (3-lb) and requires that 
the missile remain on a line (beam) joining the target and a 
control point. The target tracker is located at the control 
point; therefore, avoiding the necessity of an on-board tracker. 
Because of this, a dedicated fire control system is needed from 
launch to intercept. The range of this type of guidance is 
normally less than with the other types. A speed advantage is 
required for line-of-sight guidance since no lead angle is incor- 
porated. The main advantage of this type of guidance is the 
simple avionics required to maintain the missile in the beam. 

C. PROPORTIONAL GUIDANCE 

A proportional guidance law is one in which the rate of change 
of the missile heading is made proportional to the rate of change 
of the line-of-sight between the missile and the target. This 
is illustrated in Figure (3-lc). Since the guidance law antici- 
pates the target's future position, it can attain a higher degree 
of. responsiveness than other guidance laws. In proportional gui- 
dance the rate of change of the line-of-sight must be sensed 
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on-board the missile. Because of this requirement, and the need 
to provide anticipated steering commands, the avionics required 
are the most complex of the three guidance systems. 

D. COMPARISON OF GUIDANCE LAWS 

In early design considerations two parameters of interest 
are acceleration required of the missile and the miss distance 
attainable. Of the three guidance laws only the proportional 
law can respond to fast maneuvering targets. Since the missile 
must stay in the line-of-sight for a beam rider system, any tar- 
get maneuver will cause large excursions in the missile flight 
path, resulting in large normal accelerations. The pursuit gui- 
dance law causes similar large excursions near intercept due 
to the velocity vector always pointing at the target. 

Several system parameters affect the miss distance attain- 
able with a particular guidance law. An excellent source on 
the effects of these parameters is an article written by Dr. 
Robert Goodstein [6] . The parameters studied for their effect on 
miss distance were: 



1) 


Sensor 


Bias Angle 


2) 


Noise 




3) 


Target 


Heading 


4) 


Target 


Acceleration 


5) 


Target 


Speed 


6) 


Wind Gusts 



The results have been reproduced and are included in Figures (3-2) 
through (3-7) . Table 3-1 provides overall guidance in the selection 
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of a guidance law and is also reproduced from the above reference. 
A first glance would indicate that proportional guidance is a 
proper choice for all cases. It must be kept in mind, though, 
that cost and simplicity are also driving factors in the design 
process. Furthermore, it can be seen that, while proportional 
guidance with a high gain has good performance against maneuvering 
targets, any noise in the system will highly degrade this perfor- 
mance. For this reason another guidance law may be desired, or 
a compromise in the gain selection may have to be made in which 
some performance is given up in order to deal with a noisy system. 
A reasonable range of proportionality constants that gives good 
performance against both maneuvering targets and noisy systems 
is k = 2 to k = 6 . 

Once a guidance law is selected, a more detailed analysis has 
to be performed to determine if the maximum acceleration required 
of that particular guidance system is within the attainable 
maneuverability limits for the missile. The maximum acceleration 
required, in turn, determines the lifting surface area needed. 

A good estimate of maximum acceleration, which keeps the miss 
distance less than 50 feet, is to set it equal to three times 
the target acceleration plus ten. 

a = 3a. + 10 
m t 

Figure (3-8) shows the miss distance sensitivity to target 
acceleration. 
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TARGET HEADING SENSITIVITY 
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Figure (3-2) 



TARGET SPEED SENSITIVITY 
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Figure (3-3) 
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Figure (3-4) 



ANGLE BIAS SENSITIVITY 
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NOISE SENSITIVITY 
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WIND GUST SENSITIVITY 
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( 3 - 7 ) 



GUIDANCE LAW TRENDS FOR AIR TARGETS 
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Missile Acceleration Limit (G's) 



E. PURSUIT GUIDANCE (DETAILED ANALYSIS) 

As stated previously, a pursuit guidance law requires the 
missile velocity vector to always point at the target. For 
this reason the missile always ends up in a tail chase situation, 
with the maximum acceleration occurring at the end of the 
encounter. From this description the maximum acceleration of the 
missile can be determined. 




Reference 

Direction 



Figure (3-9). Pursuit geometry. 



From Figure (3-9) the time rate of change of the range, 
R = V T cos 6 - V M 

also 

3 = -V T sin 3/R 
or 



R 



dJ3 

dt 



V„ 



sin 



6 



R, is 
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« - A dt 

sin 3 R 



( 1 ) 



dt = 



dR 



V T cos 3 - V M 



( 2 ) 



Substituting equation (2) into equation (1) gives. 



V„ 



V, 



d3 _ T dR T 4.4. ■ , _ W M 

sin 3 ' R (V T cos 3 - V M ) ; Lettin 9 k V(] 



(cos 3 - k) AO _ dR 
sin 3 dB - R 



(3) 



Integrating equation (3) yields, 



JlnR = kJln(tan £-) - Jin (sin 3) + £nc^ 



JlnR = Jin 



c 1 tar^ 3/2' 
sin 3 



By trigonometric identity. 



(tan S/2) k = <sin B) 



(1+cos 3) 



Therefore , 



JlnR = Jin 



c 1 (sin 3) 
(1+cos 3)* 



k-1' 



and 



R = 



(sin 3) 
(1+cos 3)* 



k-1 
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From the initial condition 3 = Bq when R = R^ , 



c 



1 



R 



o 



( 1+cos 8 Q ) k 

/ * O \ 1 

(sinB Q ) 




Substituting the above equation for R into the equation for 3 
6 = -V t sinB/R 

8 _ ^t. / 1+cosB 

' R 0 \ 1+cos6 0 

The missile acceleration can be expressed as a normal and a 
tangential component, 

a = V 3 n + V t 
M M 



(sinB) 



2-k 



(sinB Q ) 



i-k 



The normal component is a m , where, 



m 



= V 3 
M p 



V V 
M T 



m 



/ 1+cos 3 
y 1+cos ^ 



) 



k 



sin 8 
sing 

0 



2-k 

1-k 



The terminal acceleration for a pursuit guidance law will occur 
at the end of the encounter (8 0). From the above expression 
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the terminal acceleration can be evaluated 



a 



0 



for 1 < k < 2 



m 




2 



2 



a 



m 



1+cosBq 



sinB^ for k = 2 



a 



00 



for k > 2 



m 



Since pursuit guidance always ends up in a tail chase situation, 
the missile will never intercept if k < 1 . Thus, for pursuit 
guidance operating against a non-maneuvering target, the velocity 
ratio should be between one and two. These results indicate 
this guidance system would not be effective against air targets; 
therefore, results for a maneuvering target were not pursued. 

F. LINE-OF-SIGHT GUIDANCE (DETAILED ANALYSIS) 






Z 



Reference 

Direction 



n 



Figure (3-10). Line-of-sight geometry 
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Figure (3-10) illustrates the geometry used to derive 
the beam rider equations of motion. The basic concept of beam 
rider guidance is that the missile is maintained in the line-of- 
sight of the target and a control point. This can be expressed 
in an equation as follows: 



<J> - 



V T sin a fc 



V„. sin a 
M m 



m 



( 1 ) 



where, r = range from point 0 to the target 
r m = range from point 0 to the missile 

From equation (1) 

r.V,, sin a = r V_ sin a. (2) 

t M m m T t 

As in the case of pursuit guidance, the missile and target acceler- 
ations can be divided into normal and tangential components. If 
the target is limited to contant g turns, and the normal component 
of missile acceleration is of interest; then. 




Differentiating equation (2) with respect to time yields, 

r.V,. sin a + r. V,, a cos a = r V_ sin a, + V_ a. cos a. 
t M mtMm mmT t T t t 



Solving for , 



a = 



m r. V„cos^ 
t M m 



— [r V m sin a. + r V_ a. cos a. - r. V M sin a ] 
a m T t m T t t t M m 
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From the original figure. 



Also 



r t = V T cos a t 



r = cos a 
m M m 



e t = <* t + * 



e t = ° t + * 



0 = a + (p 

m m 



9 = a + <J> 

m m 



The target and missile accelerations (normal components) become. 



“t = V T 



a m V M 0 m 



Collecting equations; 



e t - a t /v t 



♦ = V T sin ct t /r t 



a. = 



9 t - * 



r, = 



V T cos a fc 



r = 
m 



V cos a 
M m 



a = 



m 



= rXcostT' [*m V T sin a t 
t M m u 



+ r m V T “t cos “t 



- r, V.. sin 
t M 



° m ] 



9 = a + <b 

m m T 



a = V„ 9 
m Mm 
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The above equations are the equations of motion which describe 
the target and missile trajectories. These equations cannot be 
solved analytically except for highly specialized cases. The com- 
plete set of equations can be solved using a numerical integra- 
tion technique. If Euler's one step method is used, the algorithm 
is as follows; 



9 t (i + l) 
<t> (i+1) 

a t (i+1) 

r t (i+1) 

r (i+1) 
m 

a (i+1) 
m 

9 (i+1) 
m 



e t (i) + At 0 t (i) 

<f> (i ) + At (f>(i) 

a t (i) + At a t (i) 

r t (i) + At r t (i) 

r (i) + At r (i) 
m m 

a (i) + At a (i) 
m m 

9 (i) + At 0* (t) 
m m 



With initial conditions; 



r t (°) = r Q 



r (0) 
m 



= 0 



4,(0) = <f> 



o 



ct (0) 
m 



= 0 



e t (0) 



= 0 



9 (0) 

m 



= 9 



m. 



a t (0) 



= a. 
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The target and missile positions can be expressed as follows; 



x (i+1) 
m 


= X (i) 

m 


+ 


At 


V M 


COS 


0 

m 


(i) 


y m (i+1) 

m 


■ y m (i > 


+ 


At 


V M 


sin 


8 

m 


(i) 


x t (i+l) 


= x t (i) 


+ 


At 


V T 


cos 


CD 

r+ 


(i) 


y t (i+D 


= y t (i) 


+ 


At 


V T 


sin 


-P 

CD 


(i) 



Where 

x (0) = y (0 ) = 0 
m m 

y t (0) = r 0 cos ^0 

y t (0) = r o sin ^0 

The above equations have been programmed on the HP 9830 com- 
puter. Table 3-II is a listing of this program. The program 
asks the user for the initial conditions and the target acceler- 
ation. It also asks for the integration step increment, At. 

It should be kept in mind when using the program that the error 
involved in integrating is of order At. The output is a plot of 
missile and target trajectories as well as the missile maximum 
acceleration and time of flight. Three examples follow which 
demonstrate possible uses of the program. (Note: All angles 

are input in radians.) 
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TABLE 3- I I 



10 PRINT "THIS PROGRAM DETERMINES THE ftCCELARRTION OF" 
20 PRINT “A BEAM RIDER MISSILE AND PLOTS THE TRAJECTOR' 
30 PRINT 

40 PRINT "INPUT TIME INCREMENT FOR INTEGRATION" 

59 INPUT D 1 

60 D I M XC 250 'J , YC 250 3 ? I.JC 250 3 , VC 250 3 , AC 250 3 , RC 250 3 
70 PRINT "INPUT INITIAL TARGET RANGE" 

80 INPUT R 1 

90 PRINT "INPUT TARGET SPEED” 

100 INPUT VI 

110 PRINT "INPUT MISSILE SPEED" 

120 INPUT V 2 

130 PRINT "INPUT INITIAL LINE OF SIGHT ANGLE" 

140 INPUT PI 

150 PRINT "INPUT MISSILE ALPHA" 

160 INPUT A 2 

170 PRINT "INPUT TARGET ALPHA" 

180 INPUT A 1 

190 PRINT "INPUT MISSILE THETA" 

200 INPUT T2 

210 PRINT "INPUT TARGET THETA" 

220 INPUT T 1 

230 PRINT "INPUT TARGET ACCELERATION" 

240 INPUT G 1 

250 1=1 

260 R2=0 

265 RC 1 ]=R1 

270 AC 1 3=0 

280 XC I 3=0 

290 YC I 3=0 

300 IJC I 3=R 1 *COS < P 1 > 

310 VC I 3=R1*SINCP1> 



320 AC I 3=0 

330 PRINT " XM YM XT 

3 9 5 P PINT " D 0 V 0 U W A N T A P R I N T 0 F T H E 0 U T P U T ? O = Y E S j 1 = 

396 INPUT T9 

397 IF T9= 1 THEN 410 

400 PRINT XC I 3.i YC I 3> UC I 3? VC I 3 

410 1=1+1 

420 T 3 = G 1 ••••■ V 1 

430 P 3 = V 1 * S I N < A 1 ) ••••' R 1 

440 A3=T3-F'3 



450 R3=V1+C0S( A1 > 

460 R4=V2*C0S< A 2 > 

461 D 6 = R 4 * V 1 * S I N ( A 1 > 

4 6 2 D 7 = R 2 * V 1 * A 3 * C 0 S < A 1 > 
4 6 3 D 8 = - R 3 * V 2 * S I N ( A 2 ) 
464 D9=R 1 * V 2 * C 0 3 ( A 2 ) 
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TABLE 3-II (cont) 



465 R 4 = D 6 + Ei 7 + Ii y ,i' ■■ D 9 

466 r 4 = R 4 + P 3 
470 T 1 = T 1 + D 1 * T 3 
480 F'1=P1 + D1*P3 
490 fl 1 =fl 1 +D 1 *R3 

495 R 1 =R 1 +D 1 -*R3 

496 R2=R2+D 1 *R4 
500 R2=R2+D 1 *A4 
510 T2=T2+D 1 *T4 
520 RE I ]=V2*T4 

5 3 0 IF R B S ( R C I 1 ) < R B S ( R [ I - 1 ] > T H E H 6 0 0 
540 M= RE I ] 

600 T2=R2+P 1 
610 RE I 1=R 1 -R2 



6 2 0 X C I 1 = X [ I - 1 ] ■ + B 1 * V 2 * C 0 S ( T 2 > 

6 3 0 Y C I ] = Y E I - 1 ] + D 1 * V 2 * S I N < T 2 ) 

640 UC I ]=UC I - 1 ] + D 1 * V 1 *C0S < T 1 ) 

6 5 0 V C I ] = V C I — 1 ]+Dl*V 1 * S I N ( T 1 ) 

660 IF REID-0 THEM 397 
665 T5=Iil*<: 1-1) 

670 PRINT "INPUT MINIMUM X VALUE" 

680 INPUT X6 

690 PRINT "INPUT MRXIMUM X VALUE " 

700 INPUT X7 

710 PRINT "INPUT MINIMUM Y VALUE" 

720 INPUT Y 6 

730 PRINT "INPUT MAXIMUM Y VALUE" 

740 INPUT Y 7 

8 y 0 8 c R L E X 6 ? X 7 Y 6 j Y 7 

810 PRINT "HRS AXIS BEEN DRAWN:. 0=YES 

820 INPUT Fl 



830 IF F1=0 THEN 860 



840 XflKIS 0 j X7-" 1 0 
850 YRXIS 0> Y7/10 



8 6 0 P E H 

870 FOR W= 1 TO I 
880 PLOT XE W 1 Y E W ] 

890 NEXT W 
900 PEN 

910 FOR 8=1 TO I 
920 PLOT UESDVES3 
930 NEXT S 
940 PEN 

950 PRINT 

951 PRINT 

952 IF RE ID-RE 1-1] THEN 954 

953 M=AE I - 1 ] 

954 PRINT "THE MRXIMUM ACCELERATION I 
960 PRINT "THE TIME TO INTERCEPT IS 
1000 STOP 



1 . Example I (Non-maneuvering Crossing Target) 




Figure (3-11). Non-maneuvering crossing geometry. 



r = 359 meters 
0 

<f> 0 = 109.57° = 1.9124 rad 

0 = 17.57° = .3067 rad 

fc 0 

0 = 109.57° = 1.9124 rad 

m 0 

a = -91.0° = 01.5882 rad 
fc 0 

a t = 0 

V., = 37 3 m/sec 
M 

V T = 221 m/sec 

Table 3-III is the computer output. As indicated the missile 
maximum acceleration is, 

a = -459.25 m/sec/sec = -46.86 g's 
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TABLE 3-III 



INPUT TIME INCREMENT FOR INTEGRATION 

INPUT INITIAL TARGET RANGE 

INPUT TARGET SPEED 

INPUT MISSILE SPEED 

INPUT INITIAL LINE OF SIGHT ANGLE 

INPUT MISSILE ALPHA 

INPUT TARGET ALPHA 

INPUT MISSILE THETA 

INPUT TARGET THETA 



INPUT TARGET ACCELERATION 

I.ILJ • I..I L J 



KM 


YM 


T 

! 


YT 


DO YOU WANT A 


PRINT OF THE CM. 


J 1 P U T ? U = Y E S 1 = N U 




0 


0 


- 1 2 0 „ 2 6447 3 4 


338.256495 


-5. 1549847 IS 


17. 92341018 


-109. 7301208 


341 . 5926479 


- 9 . 197499767 


36. 1380181 


- 99. 19576 8 1 8 


344.9288007 


-12. 11 S3 1881 


54.5581 97 1 3 


— 8 8 » 66141559 


348 . 2649534 


-13. 90888892 


73.1 1385836 


-78. 127063 


351 . 601 1062 


-14.58128543 


91.75172791 


— 6 1 ' . 592 1 - ' 1 0 4 


354.937259 


-14. 14756982 


110. 3 9 6 6 8 5 9 


-57. 05835782 


358.2734118 


-12. S29971SS 


128. 9848380 


— 4 6 . 52 4 U 0 5 2 3 


3 61 . 6 0 9 5 6 4 6 


-10. 0574 0 5 8 


147.456557 


- 3 5 . 9 8 9 6 5264 


364.9457174 


— S . 4 S 4 9 9 S 4 4 8 


165.7572973 


-25.45530805 


368.2818702 


- 1 . 893028852 


183.8382156 


-14.92094745 


371.618023 


3. SI 4079320 


201 • 65 6 5 8 4 4 


— 4 • 3 8 6594 8 6 3 


374.9541758 


1 0 . U0880S33 


219. 1760013 


6. 147757728 


378 . 2903286 


17.24154285 


236 . 3664062 


1 6' . 6 8 2 1 1 U 8 2 


381 , 6264814 


25. 2S 1535 19 


253.2039305 


27.21646291 


384.9626342 


34. 01774S32 


269 . 6706019 


37.7508155 


O O Q O "7 O "7 

_■ O «»• II u_ .7 '.j i i 


48. 4 5 9 S 1 0 9 9 


285. 7539377 


48. 28516809 


39 1 . 6349398 


53. 537S7559 


30 1 . 4464555 


58.81952068 


394 „ 97 1 0926 


64. 2041 194 


316.745131 


69. 35387327 


398.3072454 


75.41315987 


331 . 6508318 


r' 9 « 3 8 8 2 2 5 8 r' 


401. 6433982 


87. 12134908 


346. 1677463 


90.42257846 


404.979551 


99.28777153 


360 . 3028325 


100. 956931 


408. 3157038 


111.8741571 


374 . 0652959 


111.4912836 


411.651 8 5 6 6 


124. 8449199 


3 8 1 7 . 4 6 6 1 U 9 7 


122. 0256362 


414. 9880094 


138. 1671393 


400.5175836 


132.5599888 


418.3241622 


151.8104937 


413.232983 


143. 0943414 


421.6 6 U 3 1 5 


INPUT MINIMUM 


X VALUE 






INPUT MAXIMUM 


X VALUE 






INPUT MINIMUM 


Y VALUE 






INPUT MAXIMUM 


Y VALUE 






HAS AXIS BEEN 


DRAWN 0=YESj 


1 = NQ 




THE MAXIMUM A 


CCELERATIOH IS 


-459.2526072 


METERS /SEC 


THE TIME TO 


INTERCEPT IS 


1 . 3 SEC 





SEC 
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450+ Target trajectory 
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Figure (3-12) . Target and missile positions for 

non-maneuvering crossing target. 



The large acceleration is typical of line-of-sight due to the 
missile requirement to stay in the beam. The trajectories are 
plotted in Figure (3-12) . 

2. Example II (Effect of V M ) 

This example is presented to study the effect of missile 
velocity on a crossing target (non-maneuvering) . 




/ 

/ 



/ 




Figure (3-13) . Crossing target geometry. 



The initial conditions are as follows; 
r Q - 4000 ra 

4> = 45° = .7854 rad 

6 = 45° = .7854 rad 

m Q 

a = 135° = 2.3562 rad 

t 0 

0 = 180° = 3.1416 rad 

t 0 

V T =200 m/sec 
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TABLE 3- IV 



THIS PROGRAM DETERMINES THE ACCELARATION OF 
A BEAM RIDER MISSILE AND PLOTS THE TRAJECTORY 

INPUT TIME INCREMENT FOR INTEGRATION 

INPUT INITIAL TARGET RANGE 

INPUT TARGET SPEED 

INPUT MISSILE SPEED = Yoo 

INPUT INITIAL LINE OF SIGHT ANGLE 

INPUT MISSILE ALPHA 

INPUT TARGET ALPHA 

INPUT MISSILE THETA 

INPUT TARGET THETA 

INPUT TARGET ACCELERATION 

XM YM XT 

D 0 Y 0 U N ANT A P R I N T 0 F THE 0 IJ T P I J T , 0 = Y ES ? 1= N 0 
INPUT MINIMUM X VALUE 
INPUT MAXIMUM X VALUE 
INPUT MINIMUM Y VALUE 
INPUT MAXIMUM Y VALUE 
HAS AXIS BEEN DRAWN? O-YES? l=NO 



THE MAXIMUM ACCELERATION IS 58-03129548 
THE TIME TO INTERCEPT IS 8 SEC 



THIS PROGRAM DETERMINES • THE ACCELARATION OF 
A BEAM RIDER MISSILE AND PLOTS THE TRAJECTORY 

INPUT TIME INCREMENT FOR INTEGRATION 
INPUT INITIAL TARGET RANGE 
INPUT TARGET SPEED 

INPUT MISSILE SPEED - ^ OO / s £- c 

INPUT INITIAL LINE OF SIGHT ANGLE 

INPUT MISSILE ALPHA 

INPUT TARGET ALPHA 

INPUT MISSILE THETA 

INPUT TARGET THETA 

INPUT TARGET ACCELERATION 

XM YM XT 

D 0 Y 0 U W ANT A P R I N T 0 F THE 0 U T F' IJ T ? 0 = Y ES j 1= N 0 
INPUT MINIMUM X VALUE 
INPUT MAXIMUM X VALUE 
INPUT MINIMUM Y VALUE 
INPUT MAXIMUM Y VALUE 
HAS AXIS BEEN DRAWN ? 0=YES> l=NO 



THE MAXIMUM ACCELERATION IS 72.78774168 
THE TIME TO INTERCEPT IS 5.6 SEC 



YT 



METERS/SEC/SEC 



YT 



METERS •■•'SEC- "SEC 
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TABLE 3- IV (cont) 



THIS PROGRAM DETERMINES THE ACCEL AR AT I ON OF 
A BEAM RIDER MISSILE AND PLOTS THE TRAJECTORY 



INPUT TIME INCREMENT FOR INTEGRATION 
INPUT INITIAL TARGET RANGE 
INPUT TARGET SPEED 

INPUT MISSILE SPEED = &OO ™/s£C 

INPUT INITIAL LINE OF SIGHT ANGLE 

INPUT MISSILE ALPHA 

INPUT TARGET ALPHA 

INPUT MISSILE THETA 

INPUT TARGET THETA 

INPUT TARGET ACCELERATION 

KM YM XT YT 

D 0 Y 0 U M A N T A PR I N T 0 F T H E 0 U T P U T , 8 = Y E S ? 1=N 0 
INPUT MINIMUM X VALUE 
INPUT MAXIMUM X VALUE 
INPUT MINIMUM Y VALUE 
INPUT MAXIMUM Y VALUE 
HAS AXIS BEEN DR ANN ? 8= YES* l=NO 



THE M A X I M U M A C C E L E R A T 1 0 N I S 8 7 . 2 2 9 5 8 8 7 1 M E T E R S / S E C / S E C 

THE TIME TO INTERCEPT IS 4.4 SEC 
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y (meters) 3000 
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Figure (3-14). Effect of V on the encounter. 



The program was run three times for a missile velocities of 
400, 600, and 800 m/sec. Table 3-IV contains the program out- 
puts and Figure (3-14) is the plot of the trajectories. As 
can be seen from the output, the missile maximum acceleration 
increases with increasing missile speed. For example, 

V/V T = 2 gives the smallest acceleration, although the maximum 
acceleration for V M / V T = 4 is not exceedingly large for this 
scenario . 

3 . Example III (Maneuvering target) 

In this example the effect of a target maneuver is inves- 
tigated. If at the time of launch the target initiates a 7 g 
(68.6 m/sec/sec) turn, the following encounter would result: 




Figure (3-15). Maneuvering target. 
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TABLE 3-V 



THIS PROGRAM DETERMINES THE ACCELARAT I ON OF 
R BERN RIDER MISSILE RHD PLOTS THE TRAJECTORY 

INPUT TIME INCREMENT FOR INTEGRATION 

INPUT INITIAL TARGET RANGE 

INPUT TARGET SPEED 

INPUT MISSILE SPEED 

INPUT INITIAL LINE OF SIGHT ANGLE 

INPUT MISSILE ALPHA 

INPUT TARGET ALPHA 

INPUT MISSILE THETA 

INPUT TARGET THETA 

INPUT TARGET ACCELERATION 

KM YM KT 

D 0 V 0 U W ANT A P R I N T OF THE 0 U T P I J T j 0 = Y ESj 1= N 0 
INPUT MINIMUM N VALUE 
INPUT MAXIMUM X VALUE 
INPUT MINIMUM Y VALUE 
INPUT MAXIMUM Y VALUE 
HAS AXIS BEEN DRAWN ? Q=YES* l=NO 



THE MAXIMUM ACCELERATION IS -355.8215593 
THE TIME TO INTERCEPT IS 5.8 SEC 



YT 



METER 
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y (meters) 10,000 




67 



Figure (3-16). Effect of a 7g maneuver upon the encounter. 



The initial conditions are; 



r 

0 


= 10000 


m 


V T = 821.436 m/sec 


*0 


= 90° = 


1.5708 rad 


V._ = 9 85.723 m/sec 
M 


0 


= 90° = 


1.5708 rad 




"b 








“t 


= -180° 






b) 








9 t 


= -90° 














a t 


= 68.60 


m/sec/sec 




From the 


output 


(Table 3-V) , 


notice the large missile acceler 


ation required (211. 


55 m/sec/sec) 


to intercept a maneuvering 



target. The trajectories are plotted in Figure (3-16). 
G. PROPORTIONAL GUIDANCE (DETAILED ANALYSIS) 




Figure (3-17). Proportional guidance geometry. 
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Proportional guidance automatically establishes a lead angle 

l 

and reacts to a changing line-of-sight . The basic guidance law 
used equates the rate of change of the missile heading to a 
constant times the rate of change of the line-of-sight. From 
the above figure this law can be expressed as. 



0 = k a 

m 

From the figure, the rate of change of the line-of-sight, a, is 
given by, 

V m sin 8 . - V„ sin 8 

• T t M m 



As in the case of pursuit and line-of-sight guidance, the para- 
meter of interest here is the normal acceleration; therefore, the 
missile and target tangential accelerations are assumed to be 
zero. In this case. 



a = 9 
m Mm 



a t = V T 0 t 



From the guidance law. 



a = V M k a 
m M 



and 



a r = V k (V sin 8 - V sin 6 ) 

m M T t M m 
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V = 0 , the time derivative of this equation is. 



Since V = 
m 



r a + r a 
m m 



k V [(V 6. cos 0 - V 0 cos 0 ) ] 

m t t t mm m 



( 1 ) 



From the original figure. 



9 = 6 + o 

m m 



m 

0 = 0 - a = 0 p 

mm m k 



0=0 (1-i) 

m m k 



Also , 

t * ► 

B t = e t - a 



» _ i 5! 

B t - 6 t k 



Making these substitutions, equation (1) becomes, 

* • * * 1 

r a = -a r + k V [ V (0 -^) cos 0 - V M 0 (l-^)cos 6 ] 

m m M T t k t Mmk m 



Since , 



a = V„ 6 
m Mm 



a t ’ V T 9 t 



r a = 
m 



-a r + k a, cos 6 , - k V.. cos & 

m - M t t Mm m 



- a m (V T cos B t ‘ V M 003 B m> 
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From Figure (3-17), 



r = V cos B - v M cos B m 
T t M m 



ra=-2ar+k V a cos 6. - k a cos B 

m m Mt t Mm m 



, k V„a, cosB. a 

1 Mt t m r « • , „ a . 

a m = r ~ [ 2r + k VjyjCosB^ ] 



m 



Collecting equations; 



r = v t cos 6 t ' V M cos B m 



8 t “ 6 t - « 



B_ = 



m 



0 - a 

m 



a = 



V T SinS t " V M sinS m 



t V, 



9 = 

m 



T 
k a 



The above equations are the equations of motion for a missile 
using proportional navigation assuming constant missile and 
target speeds. As with the line-of-sight equations, the motion 
is quite complex. The equations cannot be solved analytically 
except for' special cases. One such case will be investigated 
here. That is for a non-maneuvering crossing target. 
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1 . Example IV (Non-Maneuvering Crossing Target) 






= 0 ) 



Figure (3-18). Non-maneuvering target. 



For this case, r = -V,_ cos 8 from equation (1) 

M m 



a = -- S [2r - kr] 
m r 



r 2 - 1 <*- 2 > 

m 



Jina = (k-2) &nr + J2nc. 
m ± 



If 



a = a„ at r = r. 
m d 0 



k-2 



a = a (— ) 
m 0 r Q 



From this equation. 



if k > 2 



m 



0 as r — ► 0 
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if k 



2 



constant 



a = 
m 



if k < 2 a — * °° as r — > 0 
m ' 



From the above example it can be realized that the pro- 
portionality constant, k, must be greater than two. A more 
general analysis of the equations of motion can be obtained by 
solving the equations numerically. The same Euler's one step 
method is used here with initial conditions, at t = 0 , 



a(0) 

m 

r (0) 
B t (0) 

6 m (0) 

0 ( 0 ) 

V 0) 

9 m (0) 




0 



m 



0 



The algorithm used is as follows: 



r(i+l) = r(i) + At r(i) 

a (i+1) = a (i) + At a (i) 
mm m 

cr(i+l) = a(i) + At a(i) 

6 (i+1) = 0 (i) + At 0 (i) 
mm m 
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0. (i+1) = 6. (i) + At 0. (i) 

t t t 

8 (i+1) = 6 (i) + At 3 (i) 
mm m 

e t (i+i) = e t (i) + At 5 t (i) 



The trajectory of the missile and target can be deter- 
mined assuming the missile is at the origin at t = 0 . From the 



original figure the missile and 

x (i+1) = x (i) + At V.. 

mm M 

y (i+1) = y (i) + At V 
m u m M 

x fc (i+1) = Y t (i) + At V T 

y t (i+l) = y t (i) + At V T 

With initial conditions. 



target positions are given by, 

cos 0 (i) 
m 

sin 0 (i) 
m 

cos 0 (i) 
sin 9 fc (i) 



x (0 ) = y (0) = 0 
m m 



x t (0) 

y t (o) 



COS CT 

sin a 



0 

0 



Evaluating the missile initial acceleration can be more complicated. 

One procedure which is both realistic and of interest is to have 

the missile and target on a constant bearing - decreasing range 

course (a fc = 0) when the target initiates a constant g turn at 

t = 0. In this case a (0) =0 and the subsequent motion can be 

m 

found. 
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TABLE 3-VI 



10 PRINT "THIS PROGRAM FINDS THE MAXIMUM" 
20 PRINT "MISSILE ACCELERATION FOR A" 

30 PRINT "PROPORTIONAL NAVIGATION SYSTEM" 
35 PRINT 

3 6 D I M X [ 2 5 0 ] » Y C 2 5 0 ] , I J C 2 5 0 ] » V C 2 5 O 1 
38 D I M AC 250 1 j TC 250 ] » RC 250 3 
40 PRINT "INPUT TIME INCREMENT" 

50 INPUT D 1 

60 PRINT "INPUT NAVIGATION CONSTANT" 

70 INPUT K 



80 PRINT 


'INPUT MISSILE VELOCITY" 


90 INPUT ' 


n 




1 00 


PRINT 


"INPUT 


TARGET VELOCITY" 


110 


INPUT 


\ I 1 
1 U- 




120 


PRINT 


"INPUT 


TARGET ACCELERAT 


130 


INPUT 


A 1 




140 


PRINT 


"INPUT 


INITIAL MISSILE 


1 5 0 


INPUT 


AC 1 3 




160 


PRINT 


"INPUT 


INITIAL. RANGE" 


170 


INPUT 


RC 1 3 




180 


PRINT 


" I NPUT 


BETA TARGET" 


200 


INPUT 


B1 




210 


PRINT 


"INPUT 


BETA MISSILE" 


220 


INPUT 


B 2 




230 


PRINT 


" I NPUT 


THETA TARGET- 


240 


INPUT 


T 1 




250 


PRINT 


" INPUT 


THETA MISSILE" 


260 


INPUT 


T2 




270 


PRINT 


" INPUT 


SIGMA" 


280 


INPUT 


SI 




290 


1 = 1 






295 


PRINT 


MISSILE POSIT 


d 9 b 


PRINT 






297 


PRINT 


" 


X 1 Y 


300 


XC I 3= 


0 




3 1 0 


YC I 3= 


0 




320 


iJC I 3= 


RC 1 3* CO 


S < S 1 > 



TARGET F'OS 



'1 



330 VC I 3=R[ 1 3*SIN<S1> 

3 4 0 P R I N T "Ii 0 Y 0 U W ANT A P R I N T 0 F THE 0 U T P U T 5 O = Y E S < 1 = N 0 ' 

345 INPUT Q 1 

350 IF Q 1 = 1 THEN 365 

360 PR I NT XC I 3 j YC I 3 j U C I 3 j VC I 3 

365 1=1+1 

370 )<l 1 3 = :•••=: C 1-1 3 +V 1 *COS ( T2 > *Ii 1 
3 8 0 Y C I 3 = Y C I - 1 3 + V 1 * D 1 * S I N ( T 2 ) 

3 90 U C I 3 = U C I - 1 3 + V 2 * D 1 * C 0 S < T 1 > 

4 0 0 V C I 3 = V C I - 1 3 + D 1 * V 2 * S I N < T 1 ) 

410 R 2 = V 2 * C 0 S < B 1 > - V 1 * C 0 S < B 2 ) 

420 A 8 = K * V 1 *A 1 * C 0 S < B 1 ) / R C I - 1 3 
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TABLE 3-VI (cont) 



430 fl 9 = A C I - 1 ] * < 2 * R 2 + K * V 1 * C 0 S < B 2 > ) / RC 1-1 3 

440 R3=R8-R9 

445 TC I 3= C I - 1 > *B 1 

450 S2= < V 2 * S I H ( B 1 > - V 1 * S I H < B 2 > ) / R C I - 1 ] 

468 T'3=R 1 2 

47Ei T9=K*S2 

480 B3=T8-S2 

490 B4=T9-S2 

500 R[ I ]=R[ 1-1 3+Bl*R2 

510 R[ I ] = RC 1-1 ] + D 1 *A3 

520 S1=S1+B1*S2 

530 T 1 =T 1 +B 1 *T3 

540 T2=T2+B1*T9 

550 B1=B1+B1*B3 

560 B2=B2+B 1 *B4 

570 IF RCIKRCI-1] THEN 600 

580 Z 1 =AC I ] 



600 


IF RC I ]>0 THEN 358 






601 


IF RBS < RC I 3><RBS(R[ I- 


-13) THEN 603 




602 


Z 1 =R[ 1-1 1 






603 


Z2=TC I ] 






6 1 3 


PRINT 






619 


PRINT 






620 


PRINT "MAXIMUM MISSILE ROC IS"Z1" 


METERS/SEC 


640 


PRINT "MISSILE TIME i 


OF FLIGHT IS 


" Z2 " SEC" 


700 


PRINT "INPUT MINIMUM 


VALUE OF X" 




710 


INPUT X6 






720 


PRINT "INPUT MAXIMUM 


VALUE OF X" 




730 


INPUT X7 






740 


PRINT "INPUT MINIMUM 


VALUE OF Y" 




750 


INPUT Y6 






760 


PRINT "INPUT MAXIMUM 


VALUE OF Y" 




770 


INPUT Y7 






730 


S C ALE X 6 j X 7 > Y 6 > Y 7 






790 


PRINT "HAVE THE AXIS 


BEEN DRAWN) 


0-YES ? 1 =NO " 


791 


INPUT Fl 






792 


IF F 1 =0 THEN 838 






800 


XRXIS 0 j X 7 1 0 






3 1 0 


YAXIS 0i. Y 7/ 10 






330 


PEN 






340 


FOR J=1 TO I 






350 


PLOT XU J ] ? YC J ] 






360 


NEXT J 






870 


PEN 






380 


FOR M=1 TO I 






890 


PLOT UCW3)VCW3 






900 


NEXT W 






910 


PEN 







1000 STOP 
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The above algorithm was programmed on the HP 9830 computer. 
The listing is included in Table 3-VI. The inputs and output 
of the program are the same as the line-of-sight guidance pro- 
gram. Several examples follow which demonstrates the use of the 
program. 

2. Example V (Crossing Maneuvering Target) 




\ 

\ 

\ 




V M = 


373 m/sec 


CD 

ft 

O 

II 


17.57° = 0.3067 


rad 


II 

> 


221 m/sec 




-91.0° = -1.5882 


rad 


a t 


156.8 m/sec/sec 




-38.0° = -0.6632 


rad 


r o 


359 m 


a = 
o 


109.57° = 1.9124 


rad 


0 


71.57° = 1.2491 rad 









m 0 

The output is listed in Table 3-VII. An important aspect of this 
problem is the maximum acceleration required • (202 . 92 m/sec/sec). 



TABLE 3-VII 



THIS PROGRAM FINDS THE MAXIMUM 
MISSILE ACCELERATION FOR A 
P R 0 P 0 R T 1 0 N A L N A V I G A T 1 0 N S Y S T E M 

INPUT TIME INCREMENT 

I H P U T N A V I G A T 1 0 N C 0 N S T A N T 

INPUT MISSILE VELOCITY 

INPUT TARGET VELOCITY 

I N P U T TARGET A C C E L E R A T 1 0 N 

INPUT INITIAL MISSILE ACCELERATION 

INPUT INITIAL RANGE 

INPUT BETA TARGET 

INPUT BETA MISSILE 



NF'UT THETA TAR 


GET 






NPUT THETA MIS 
NPUT SIGMA 


SILE 






MISSILE POS 


IT 


TARGET PO 


S I T 


X 1 


Y 1 


..... 


Y b 


1 0 Y 0 U W ANT A P R I N T 0 F T H E 0 U T P U T > 0 = Y E S ? 1 = N 0 




0 


0 


-120. 2644734 


338,256495 


5 . S 9 6 S 8 8 3 8 2 


1 7 ■ b 9 b \d 6 0 


- 1 09 „ 730 1 288 


34 1 , 5926479 


11.72925942 


b b » 4 0 7 7 b fo b 6 


-99,32072172 


345 , 3003306 


1 7 . 50036490 


53. 14238894 


-89. 04937488 


3 4 9 . 3 i 4 o r* r* r 


23.21064114 


70. 896691 10 


-78.92900524 


353.8111620 


2 y . y 5 y 6 u 4 7 s 


88 . 6709 1421 


- 6 8 . 9 7 2 34776 


358 . 60360 1 


34. 44056825 


1 y 6 . 4 6 5 9 7 4 9 


-59. 19193143 


3 6 3 . 746 1 6 4 2 


39 . 95057982 


124.2834458 


“49 . 60006345 


369 . 2323805 


45. 38038997 


142, 1255214 


-40. 20881377 


375, 0553462 


50. 71944940 


1 5 y ■ 9 4 44f-,. -! 


-3 1 , 02999988 


38 1 . 2077340 


55.95494186 


177.8950225 


-22.07517195 


387 . 68 1 8020 


61 . U 7 1 y 4 9 6 y 


195.82934 


-13.3555983 


394 . 4694036 


66. 053O71 1 


213. 8018172 


-4.882251242 


4 y 1 , b 6 1997 6 


70.87956373 


231.8164631 


3 , 3 3 4 2 0 6794 


408 . 9506590 


75.53054295 


249.8772187 


11.2 8 3 4 3 y 6 


416. 6260903 


79. 9 8 3 7 1 8 9 2 


267 , 9877605 


18.95543527 


424 . 578633 


84.21557618 


286. 1512922 


26. 3 4 O 5 4 8 { 2 


432 . f* 9 y 2 y y 2 


88.20169238 


304.370331 


33 . 42948387 


44 1 . 27 4688b 


91 . 91709127 


322, 6464989 


40. 21332036 


449 . 997 1 9 1 8 


95. 33662288 


340. 9803284 


4 y . 6 8 3 5 2 1 r 4 


458.9548139 


9 y . 4 3 5 3 6 1 3 2 


359.3710956 


5 2 . 8 3 1946 2 


468. 1 3 6 2 8 3 


101 . 1890086 


377 . 8 1 66893 


5 8 . 6 5 0 8 5 6 8 9 


477. 5308458 


103. 5742899 


396 . 3 1 3525 1 


64, 13293155 


487. 1242815 


105.5693244 


414.8565113 


6 9 . 2 71271 8 2 


496. 906917 


107. 1539546 


433 . 4390689 


74.05941 184 


506 . 8656426 


108.3100167 


452 . 0532039 


7 8 , 4 9 1 3 2646 


516.9 8 7 9 2 6 6 


109.0215346 


470. 6896263 


82.56143878 


527.2610316 


109. 2748269 


489. 3379062 


86. 26462715 


537,6720304 


1 09 . 0585 1 42 


507. 9866517 


89 . 59623 167 


548.2078223 


108,3634271 


526. 62 3 6 9 4 3 


92.55206 


558.8551496 


1 U r' , 1 y 2 4 2 y u 


545 . 2362632 


9b. 12839268 


5 y 9 . 6 y y 6142 


105. 5101054 


563 . 8 1 1 1 35 


9 7 . 3 2198776 


580 . 4306946 


1 y 3 . 3 4254 1 y 


582 . 334746 1 


99. 13008494 


591 . 3317626 


100. 6769169 


600. 7932657 


1 00. 550409 


602, 2901010 
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TABLE 3-VII (cont) 



MAXIMUM MISSILE ACC IS 20 
MISSILE TIME OF FLIGHT IS 
INPUT MINIMUM VALUE OF X 
INPUT MAXIMUM VALUE OF X 
INPUT MINIMUM VALUE OF Y 
INPUT MAXIMUM VALUE OF Y 
HAVE THE AXIS BEEN DRAWN j 



9 2 1 2 3 8 9 
" I . 7 " " " SE 



Q=YES !' 1 - N 0 



METER 
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y (meters) 1000 




i 
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Figure (3-20) . Trajectories of crossing, maneuvering intercept. 



The same problem was run for line-of-sight guidance with no 
target acceleration in Example I. The maximum acceleration 
was 459.25 m/sec/sec. This points out the advantage of pro- 
portional navigation over line-of-sight guidance. Figure (3-20) 
is a plot of the trajectories. 

3 . Example VI (Effect of k) 

This example demonstrates the effect of varying the pro- 
portionality constant, k. The scenario is as follows. 




/ 



/ 




Figure (3-21). Initial geometry, 



V = 208m/sec 



V = 413 m/sec 



a fc = 68.60 m/sec/sec 

8. = -106° = -1.85 rad 

t 0 

8 = -29° = -.5061 rad 

m Q 
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TABLE 3-VIII 



INPUT TIME INCREMENT 
INPUT NAVIGATION CONSTANT = 3 
INPUT MISSILE VELOCITY 
INPUT TARGET VELOCITY 
INPUT TARGET ACCELERATION 
INPUT INITIAL MISSILE ACCELERATION 
INPUT INITIAL RANGE 
INPUT BETA TARGET 
INPUT BETA MISSILE 
INPUT THETA TARGET 
INPUT THETA MISSILE 
INPUT SIGMA 

MISSILE POSIT 

X 1 

DO YOU WANT A PRINT 



MAXIMUM MISSILE ACC 
MISSILE TIME OF FLIGHT IS 3.2 SEC 
INPUT MINIMUM VALUE OF X 

INPUT MAXIMUM VALUE OF X 

INPUT MINIMUM VALUE OF Y 

INPUT MAXIMUM VALUE OF Y 

H A V E T H E A X I S BEE N D R A W N - 0 = Y E S , 1 = N 0 



THIS PROGRAM FINDS THE MAXIMUM 
MISSILE ACCELERATION FOR A 
P R 0 P 0 R T 1 0 N A L N A V I G A T 1 0 N S Y S T E M 

INPUT TIME INCREMENT 

INPUT NAVIGATION CONSTANT = V 

INPUT MISSILE VELOCITY 

INPUT TARGET VELOCITY 

INPUT TARGET ACCELERATION 

INPUT INITIAL MISSILE ACCELERATION 

INPUT INITIAL RANGE 

INPUT BETA TARGET 

INPUT BETA MISSILE 

INPUT THETA TARGET 

INPUT THETA MISSILE 

INPUT SIGMA 

MISSILE POSIT TARGET POSIT 

X 1 y 1 X 2 

D 0 Y 0 U W ANT A PRINT OF THE 0 U T P IJ T ? O = Y E S ? 1 = N 0 



M A X I M IJ M M I S S I L E A C C IS 71-11 3 3 3 1 9 7 METE R S / S E C 

MISSILE TIME OF FLIGHT IS 3.2 SEC 



TARGET POSIT 

Y 1 X2 

0 F THE 0 U T P U T > O = Y ES» 1= N 0 

I S 8 3.70521471 M E T E R S / S E C 
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TABLE 3-VIII (cont) 



THIS PROGRAM FINDS THE MAXIMUM 
MISSILE ACCELERATION FOR H 
P R 0 P 0 R T 1 0 N A L H A V I G A T 1 0 H S Y S T E M 

INPUT TIME INCREMENT 
I N P IJ T N A V I G A T 1 0 H C 0 N S T A N T - S’ 

INPUT MISSILE VELOCITY 

INPUT TARGET VELOCITY 

I N P U T TARGET A C C E L E R A T 1 0 N 

INPUT INITIAL MISSILE ACCELERATION 

INPUT INITIAL RANGE 

INPUT BETA TARGET 

INPUT BETA MISSILE 

INPUT THETA TARGET 

INPUT THETA MISSILE 

INPUT SIGMA 

MISSILE POSIT TARGET POSIT 

X 1 Y 1 X 2 

D 0 Y 0 U W A N T A PR I N T 0 F T H E 0 U T P U T > 0 = Y E S j 1 = N 0 



M A X I M U M M I S S I L E A C C I S 6 5 . 3 8 57724 8 METE R S S E C / S E C 

MISSILE TIME OF FLIGHT IS 3.2 SEC 

INPUT MINIMUM VALUE OF X 

INPUT MAXIMUM VALUE OF X 

INPUT MINIMUM VALUE OF Y 

INPUT MAXIMUM VALUE OF Y 

HAVE THE A X I S BEE N D R A W N > 0 = Y E S , 1 = N 0 
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y (meters) 1200 r 
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Figure (3-22) . Effect of varying 



.8483 4ad 



6 = 48.6° = 

m 0 

0 = -28.4° = -.4956 rad 

t 0 

a Q = 77.6° = 1.3544 rad 

The problem was run for k = 3,4, and 5. 

Table 3-VIII is the output. It can be seen that the 
effect of increasing the proportionality constant is to decrease 
the maximum acceleration required. Figure (3-22) is a plot of 
the trajectories. 

H. DESIGN EXAMPLE (GUIDANCE LAW SELECTION) 

From the examples given in this chapter, it can be seen that 
a missile designed to encounter a highly maneuverable target, 
such as a fighter, requires a proportional guidance law to limit 
the maximum acceleration required of the missile. To select 
the proportionality constant, it was assumed that the threat 
could maintain a constant 7g turn at = 1.5 and an altitude 
of 10,000 feet. 

Three cases were investigated, (1) A head-on encounter with 
the target initiating a turn at 10,000 meters range, (2) A 
crossing encounter in which the target turns into the missile at 
10,000 meters range, and (3) An oblique, closing encounter in 
which the target turns into the missile. The scenario and 
computer outputs are shown in Figures (3-23) , (3-24) , and (3-25) . 

From this analysis the crossing encounter requires the largest 
acceleration (263.2 m/sec/sec). 

For the crossing scenario then, the missile speed was varied 
and the results indicated that as the speed increased the maximum 
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V = 821 m/sec (M = 2.5) 
m 

V = 492 m/sec (M = 1.5) 
r Q = 10,000 m 



0 = 180° 
0 



THIS PROGRAM FINDS THE MAXIMUM 
MISSILE ACCELERATION FOR A 
P R 0 F' 0 R T 1 0 N A L N A V I G A T 1 0 N S V S T E M 

INPUT TIME INCREMENT 

I N F' U T N A V I G A T 1 0 N C 0 N S T A N T 

INPUT MISSILE VELOCITY 

INPUT TARGET VELOCITY 

INPUT TARGET ACCELERATION 

INPUT INITIAL MISSILE ACCELERATION 

INPUT INITIAL RANGE 

INPUT BETA TARGET 

INPUT BETA MISSILE 

INPUT THETA TARGET 

INPUT THETA MISSILE 

INPUT SIGMA 

MISSILE POSIT TARGET POSIT 

XI Y 1 X2 

D 0 Y 0 IJ H ANT A P R I H T 0 F T H E 0 IJ T P IJ T , 0 = Y ES ? 1 = N 0 



6 =0 



m. 



6^ = 180 



a o " 

k = 



0 

4 



M A X I M I J M M I S S I L E A C C IS 126.1911 9 94 METE R S / 

MISSILE TIME OF FLIGHT IS 9 SEC 

INPUT MINIMUM VALUE OF X 



Figure (3-23). Head-on scenario. 
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CO 



7 g turn 




Tq = 10,000 m 

= 492 m/sec 

V = 821 m/sec 
m 

B. = -120° 

0 

B = -46.2° 
m_ 



9 



m c 

7 o 

k 



= 73. 8 C 
= 120 ° 
= 4 



THIo HRLi. RHT- 1 FINDS THE MAXIMUM 
MISSILE HCCELERAT I OH FUR fl 
F' R U F' U R T I Ci H fl L H fl V I G fl T I ij N ' S Y S T E M 

INPUT TIME INCREMENT 
I N F' U I N fl V I G fl T 1 0 N C fl N S T fl N T 
INPUT MISSILE VELOCITY 
I N P U I T fl R GET V E L 0 C I T Y 
I N P I J T TARGET fl C C E L E P fl T 1 0 N 

SSm initial tmU LE fl6cELERfiTI0N 

INPUT BETA TARGET 
INPUT BETA MISSILE 
INPUT THETA TARGET 
INPUT THETA MISSILE 
INPUT SIGMA 

MISSILE POSIT 



DO 



X 1 



Y 1 



target posit 



N ANT fl P R I N T OF THE 0 U T P IJ T 0 = Y E S t 1 = N 0 



T 



MAXIMUM MISSILE ACC IS 2 AS l MPTP . 

MISSILE TIME OF FLIGHT IS o * o ' l -. Pl , METEko " •- ,El -"'SEl 
INPUT MINIMUM VALUE OF x" oEu 



Figure (3-24) . Crossing scenario. 
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__ 



v t 

V m 



6 



m. 



r 



0 



492 m/sec 

821 m/sec 
135° 

53.9° 

10,000 m 




135° 

53.9° 

0 

4 



THIS PROGRAM FINDS THE MAXIMUM 
MISSILE ACCELERATION FOR A 
P R 0 P 0 R T 1 0 N A L N A V I G A T 1 0 N S V S T E M 

INPUT TIME INCREMENT 

I N P U T N A V I G A T 1 0 N C 0 N S T A N T 

INPUT MISSILE VELOCITY 

INPUT TARGET VELOCITY 

INPUT TARGET ACCELERATION 

INPUT INITIAL MISSILE ACCELERATION 

INPUT INITIAL RANGE 

INPUT BETA TARGET 

INPUT BETA MISSILE 

INPUT THETA TARGET 

INPUT THETA MISSILE 

INPUT SIGMA 

MISSILE POSIT TARGET POSIT 

XI • Y 1 X2 

D 0 Y 0 U H ANT A P R I N T 0 F THE 0 U T PUT j 0 = Y ES ? 1 = N 0 



MAXIMUM MISSILE ACC IS- 156. 9795962 METERS/SEC/SEC 

MISSILE TIME OF FLIGHT IS 8.8 SEC 
INPUT MINIMUM VALUE OF X 



Figure (3—25) . Oblique scenario. 
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TABLE 3 -IX 



INPUT TINE INCREMENT 
I N P U T N ft V I GRT 1 0 N C 0 N S T fi N T 

I H P IJ T M I S S I L E V E L 0 C I T Y ** = 2 . O 

INPUT TARGET VELOCITY 

INPUT TARGET ACCELERATION 

INPUT INITIAL MISSILE ACCELERATION 

INPUT INITIAL RANGE 

INPUT BETA TARGET 

INPUT BETA MISSILE 

INPUT THETA TARGET 

INPUT THETA MISSILE 

INPUT SIGMA 

MISSILE POSIT TARGET POSIT 

XI VI X2 

Ii 0 Y 0 U 14 A N T A P R I N T OF T H E 0 U T P U T » 0 = V E S » 1 = N 0 



MAXIMUM MISSILE ACC IS 275.5734145 METERS/SEC- 

MISSILE TIME OF FLIGHT IS 10.3 SEC 

INPUT MINIMUM VALUE OF X 



THIS PROGRAM FINDS THE MAXIMUM 
MISSILE ACCELERATION FOR A 
P R 0 P 0 R T 1 0 N A L N A V I G A T 1 0 N S Y S T E M 

INPUT TIME INCREMENT 

INPUT NAVIGATION CONSTANT _ _ 

INPUT MISSILE VELOCITY > M r-\~ ‘ 3 -° 

INPUT TARGET VELOCITY 

I N P U T TARGET A C C E L F. R A T 1 0 N 

INPUT INITIAL MISSILE ACCELERATION 

INPUT INITIAL RANGE 

INPUT BETA TARGET 

INPUT BETA MISSILE 

INPUT THETA TARGET 

INPUT THETA MISSILE 

INPUT SIGMA 

MISSILE POSIT TARGET POSIT 

V \ y j. ■’"» 

D 0 Y 0 U ii ANT A P R I N T 0 F T H E 0 U T P U T , 0 = Y ES * 1= N 0 



MAXIMUM MISSILE ACC IS 248,233583 METERS/SEC 

MISSILE TIME OF FLIGHT IS 7.8 SEC 
INPUT MINIMUM VALUE OF X 
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NRVIGRTION CONSTANT 



acceleration decreased. The results for M, 



M 



2 . 0 and = 3.0 



are shown in Table 3-IX. This result indicates the desirablity 
of retaining the missile Mach number originally selected. 



proportionality constant was then varied from k = 2 to k = 6 . 

The results are plotted in Figure (3-26). If the maximum 

acceleration is limited to 31 g's (3a fc + 10), the required pro- 

portionality constant is k = 3.75. This is well within the 
desirable range of 2 - 6 indicated earlier. 

From this analysis the required performance objectives are 



For the crossing case and a missile speed of = 3.0, the 



M 



3.0 



m 



k = 3.75 




max 



31 g's 
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IV. SIZING THE DIAMETER 



The missile diameter is determined by one of three driving 
factors. For relatively short range missiles the diameter will 
be fixed by either the warhead or the seeker requirements. As 
might be expected, for longer range missiles the diameter will 
more likely be fixed by either the warhead or the seeker require- 
ments. As might be expected, for longer range missiles the 
diameter will more likely be fixed by the propulsion require- 
ments in order to prevent excessive propulsion system lengths. 

An initial estimate of the missile diameter must be made at 
this point in order to proceed with the design. The initial 
seeker requirement can be determined from a knowledge of the 
lock-on range requirement found in Chapter 2. The warhead 
necessary to inflict a "kill" can also be estimated from infor- 
mation about the target and characteristic explosives. The 
propulsion requirement cannot be determined because of the lack 
of any aerodynamic drag or weight information at this point. For 
this reason the missile diameter will be now sized for seeker 
or warhead requirements. The missile propulsion requirements 
will be determined later in the design process, and it may be 
necessary at that time to resize the missile to meet these pro- 
pulsion requirements. 

Selection of the type of seeker depends upon the operational 
arena of the missile. The seeker of a shoulder fired, battle- 
field missile would not be the optimum seeker of a shipboard 
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missile where antenna and component sizes are not limiting 
factors. All types of guidance use some portion of the electro- 
magnetic spectrum. The three primary areas of use are the 
electro-optical, infrared, and radio frequencies. The milli- 
meter wave section of the spectrum is also of current interest 
in the design of missiles due to small component size and will 
also be discussed. The following table lists some of the major 
advantages and disadvantages of the three. 



Optical 



Infrared 



RF 



Advantages 



Target resolution (de- 
tail) 

Real time information 
Three dimension effect 

Improved resolution 
over RF 

Longest range 

Least absorption and 

attenuation 



Disadvantages 

Bad weather degrades 
Night use degrades 

Attenuation due to 
aerosols and atmosphere 

Larger components 



A. THE RADAR RANGE EQUATION 

An omnidirectional antenna is one that radiates power in all 
directions equally. If the power radiated by an antenna is P , 
the power density at a distance R fc from the source is given by, 

2 

Power density = P t /(4ir R fc ) 

2 

4tt r^ = area of a sphere of radius R fc 

Since antennas are normally directive instead of omnidirectional, 
most of the power is radiated in a particular direction. The 
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gain, G fc , is a measure of the increased power from a directive 
antenna as compared to an omnidirectional antenna. Therefore the 
power density from a directive antenna can be expressed as. 



Power density = 



P G 
t t 



4 7T R , 



This is the power density which arrives at the target. The tar- 
get intercepts a portion of this energy and reradiates it in the 
opposite direction. The radar cross section, a, is a measure of 
the effective area of the target. The power radiated by the 
target is p ec j 10 ^ 



echo 



P t G t a 
2 

4 it R. 



This energy propagates as if it were radiated by an omni- 
directional antenna. Therefore if the receiving antenna is a 
distance, R^, away the power density at the receiver is 



(Power density) 



p g a 

t t 



( 4 TT R fc 2 ) (4ir R r 2 ) 



If the energy is intercepted by the receiving antenna, which has 
an effective area as seen by the returned energy of ; then the 
power received by the radar , P , is 



94 



( 1 ) 



p = 
r 



P, G a A 
t t r 

( 4 TT R fc 2 ) (4 tt R r 2 ) 



This is the simplest form of the radar equation and can be used 
to determine the size of antenna required. 

B. ACTIVE RADAR HOMING 

Active homing is the method of missile guidance in which the 
radar transmitter and receiver are located on-board the missile. 
In this case the same antenna is used for both transmitting and 
receiving. The radar equation then becomes. 



P 

r 



P. G. 0 A. 
t t t 

2 2 
( 4 TT R^V 



where, R. = R = R and A = A, . 

t r r t 

The minimum power for which the target can be detected, P . , is 

min 

a function of many variables. A full development of this term 
can be found in reference (8) . 



P . = k T n B F ( rr— ■) 

min 0 n n N' m . 

0 min 



Boltzmans constant is k = 1.38 x 10 23 joule/°k. The value of 



kTg at room temperature is 4 x 10 



-21 



watt/cps of bandwidth. The 



bandwidth, B^, noise figure, F^, and minimum signal to noise 

ratio, (Sq/N ) are all functions of the receiver. Typical 
min 
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values are listed below. 



B = 1 MHz 
n 

F n = 7.5 db = 2.37 (for crystal mixer) 

(Sq/N.) = 14.7 db = 5.43 (for probability of detection, 

min 

Pp = 0.9 and probability of false alarm, 

Pf a = 1/15 minutes) 

-14 

The above values give P . = 5.15 x 10 watts. This is the 

mm 

value which will be used throughout this section. 

From antenna theory the gain is related to the effective 
antenna area by, 



G = 



4 TT A 



The maximum radar range can then be shown to be 



max 



V4 



P t * A t 

4 7T X 2 P . 

mm 



( 1 ) 



Equation (1) for radar range does not include any system losses. 
It also does not include the statistical nature of several of 
the parameters. Because of these assumptions the actual range 
of a radar may be as small as one^ialf of what the radar range 
equation predicts for laboratory conditions. For this reason, 
twice the required range should be used when using the above 
equation. 
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1 . Example 



The AN/APQ-153 is the airborne attack radar system used 
on the F-5E aircraft. The following parameters apply to this 
radar. 



f Q = 8-10 GHz 

g 

, 3 x 10 m/sec 

A = c/f = q — = .0333 m 

u 9 x 10 /sec 

c = speed of electromagnetic propagation 

g 

= 3 x 10 m/sec 
P = 80 kW 

Antenna = Parabolic dish 30.5 x 40.6 cm 
A t = 0.12383 m 2 

for a target of 1 square meter of radar cross-section. 



R = 36 . 16 km = 19 . 5 nmiles 
max 

For the case of an active homing radar, the size of the 
transmitting and receiving antenna is the parameter of interest. 
The antenna size may very well drive the design diameter of the 
missile. The antenna diameter can be expressed as follows. 



Let 



A. = 



C = 



2 4 

4ir X P . R 

mm max 



P t G 



4tt X P . 

mxn 



1/2 



P t a 



1/2 



(2) 



then 



rr d t 2 

A = - 4 - = C R 

t 4 max 



4C R 



max 

7T 



therefore, d = 2 R 




max\/ tt 



25 1 - 




(S3H0NI ) d313Wyia yNN31NU 



• 

Q) 

N 

•H 

fd 



a; 

p 

•H 
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MAX RRNGE (NMILES) 



Equation (2) has been plotted in Figure (4-1) for various 
values of the transmitter frequency. From this plot and a know- 
ledge of the maximum lock-on range required, the antenna size 
can be determined. 

From equation (2) and Figure (4-1) there are two obvious 
ways to decrease the antenna size required. (1) Increasing fre- 
quency is the best way to reduce antenna and electronic component 
sizes. The current trend is toward higher frequencies. (Milli- 
meter waves.) One problem is that the equations developed in 
this section do not include atmospheric attenuation. For fre- 
quencies above about 30 GHz the absoprtion due to atmospheric 
gases increases. This is shown on Figure (4-2). As indicated 
on this figure there are "windows" where the attenuation is less. 
These "windows" occur at frequencies of 34 GHz, 94 GHz, 140 GHz, 
and 220 GHz. These are the frequencies where most of the cur- 
rent research and development is going on. As the frequency in- 
creases, the wavelength approaches the size of rain droplets. 

For this reason, radar performance is greatly reduced in inclement 
weather. (2) Increasing transmitter power will also decrease 
the size of antenna necessary. The limiting factor in this area 
is the lack of high power sources. In the millimeter range the 
available power from current traveling wave tubes is 50-100 
watts. Increasing power is obviously confined to size and weight 
limitations of missile components. 
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Figure (4-2). Atmospheric Absorption [8]. 
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C. SEMI-ACTIVE HOMING 

The advantage of semi-active homing is obvious when the radar 
range equation is investigated. From equation (1) , 



2 2 
R. R 
t r 



P. G. o A 
t t r 

(4tt) 2 P . 

min 



In the above equation the missile range from the target is R^. 
The transmitting and receiving antennas are at different ranges 
and have different characteristics in this case. As before, 

4tt A fc 



therefore , 





a 



4tt X 2 



P . 
mm 



The main advantage is in the transmitter characteristics. Since 
the transmitter is not located in the missile, it is not normally 
limited in size and weight requirements. In the above equation 
if a transmitter power and standoff range, R , is chosen the 
receiving antenna can be sized for a maximum homing range, R , 
of the missile. 

1 . Example 

P t = 100 km 

R t = 100 nmiles 

f Q = 10 GHz 
2 

A fc = 4 m 

P . = 5.15 x 10 -14 W 

min 
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R „ (nmiles) d (in) 

max r 



10 


5.8 


20 


11.6 


50 


29.1 



It can be seen from comparing these numbers to those of Figure 
(4-1) that the required antenna size is less than one half of 
that required for an active homing radar of the same frequency. 

D. DESIGN EXAMPLE (ANTENNA SIZING) 

An active radar was assumed in Chapter 2 to decrease the 
required missile range. The lock-on range was 10,000 m or 
5.4 nmiles. As stated earlier, twice this number should be used 
for determining antenna size. From Figure (1) for a range of 
10.8 nmiles, and a transmitter power of 10 kW at f^ = 20 GHz, 
the required antenna size is d = 10 inches. 



E. INFRARED SEEKERS 



In the design of missile seekers two parameters of primary 
importance are range and size. The idealized range for an infra- 
red tracker relates these two factors. The idealized range is 
the range at which the signal-to-noise ratio is unity and is 
given by, 



R 0 = 



D*T T D 2 J "I 1,/2 

a IR a 



4 \/A f A 



( 1 ) 
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Typical value 



Specific detectivity 

Transmission through 
atmosphere 

Transmission through 
IR optics 

Aperture diameter 

Radiant intensity 

Receiver bandwidth 

Sensitive area of 
detector 



D* 



10 



10 



T 0-1.0 

a 

t ir °- 1 - 0 

D 

a 

J 10 3 

Af 10 3 

10 ^ - 10 ^ cm 3 



The derivation of equation (1) and its use are the subjects of 
this section. Some references (9) may give the above equation 
in terms of the Noise Equivalent Intensity, NEI. 



= UJ 



R o [nei j 



1/2 



where 



NEI = 



/Af A. 



D* A T TIR 
a a 



1. Planck's Law 

The radiant emittance of a body is a measure of the radiant 
power per unit area emitted from the surface. 

W = — watt/cm 3 (2) 



The spectral radiant emittance is the radiant emittance per unit 
wavelength interval, 
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w. 



3W . . , 2 

■gy watt/cm y 



-6 



y = micron = 10 meters. 



Planck's law gives the blackbody spectral radiant emittance as 
a function of wavelength and temperature, 



( V BB 



2 i h c‘ 



exp (hc/XkT) -1 



(3) 



-34 

Planck's constant = 6.6238 x 10 Joule-sec 



h 
c 
X 
k 

T = Absolute temperature, '"'K 



g 

speed of light = 3 x 10 m/sec 
wavelength 

Boltzmann's constant = 1.38 x 10 - ^ Joule/°K 



Figure (4-3) shows equation (3) for various absolute tempera- 
tures. As can be seen the wavelength at which maximum radiant 
emittance occurs varies with temperature. This maximum occurs 
at a wavelength given by Wien's displacement law, X 



max 



X T = 2897. 8y K 
max M 



2 . Emissivity 

Actual bodies do not emit radiation according to Planck's 
law. A more typical plot of radiant emittance is shown in 
Figure (4-4). Spectral emissivity, e^, is defined as the ratio 
of the actual spectral radiant emittance to the blackbody spectral 
radiant emittance, 

W, 



£ X (W,) 



BB 
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Spectral radiant emittance (W 




Figure (4-3) . Spectral radiant emittance. 
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A 



Figure (4-4). Actual spectral radiant eraittance. 

As shown in the above figure, e^, may vary with wavelength. A 
grey body is defined as one which has constant spectral emissivity, 

e = = constant 

IR systems normally use filters to limit the accepted 
radiation to a specific wavelength band. The radiant emittance 
of a body between wavelengths A^ and A^ becomes, 



In Figure (4-5) the surface at the orgin emits a total 
energy WA into a hemisphere normal to A. 

The radiance is defined as the radiant power per unit 
solid angle per unit projected area. 




1 



Radiance = N 



cos6 9A9ft 



_1 3 9P 

cos0 9ft 3A 
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z 




Figure (4-5) . Radiant emittance of a body. 

From equation (2) 

N = - t watts/steradian cm^ 

cos0 8ft 

The radiant intensity is defined as the radiant power per unit 
solid angle from a point source. 

3P 

Radiant intensity = J = watt/steradian 

The following is a summary of the definitions of radiant energy 
quantities , 

W = radiant emittance, watts/cm^ 

J = radiant intensity, watts/steradian 

2 

N = radiance, watts/s teradian-cm 
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From the definitions the relationship between W, J and N are. 



w = TTN = TT J/A 

From the above definitions and Figure (4-5) the energy into solid 
angle, ft^ is, 

W 

NA ft, cos 0 = — ft, A cos 9 
1 it 1 

The energy into solid angle is given by, 

W 

NAft„ = £ ft- A 
2 TT 2 



3 . Energy into a Hemisphere 

From the definition of radiance and the radiant emittance, 
the radiance in terms of the radiant emittance can be found, 




Figure (4-6). Energy emitted into a hemisphere. 
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From Figure (4-6) the incremental solid angle is 



dft 



r sin 9 d <j) r d0 

2 

r 



sin 9 d 0 d <j> 



If the above surface, A fc , is considered a Lambertian surface, 
the radiance, N, is independent of the direction of radiation. 

dW = cos 9 N dfi 



dW = cos 0 N sin 0 d 9 d <j> 



The total radiant emittance into a hemisphere above the surface 
is then, 

2tt -tt/2 



W 



■ / / N cos 0 sin 0 d 0 d <j> 



o o 



W = 2tt N 



1 ■ 2 .1 

tt sin 0 

Jo 



tt/2 



W = 7T N 
4 . Targets 

Infrared targets include a wide variety of radiation 
sources. The radiance of most bodies can be divided into that 
due to self-emission and that due to reflection of incident 
radiation. 



N = N + N 
e r 

The relative magnitude of these contributions depends on a number 
of factors and varies from target to target and operating 



environment. 
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a. Self -emission . 



Self-emission, also referred to as thermal emission, 
depends primarily on the temperature of the body and the emis- 
sivity. The most often used appraoch is to consider the body as 
a grey body which emits radiation according to the Stefan- 
Boltzmann law. 

4 

N = - — watts/steradian-cm^ 
e 7T 

where a = 5.67 x 10 12 watts/cm 2 (°K) ** 

This term is the total radiance (over all wavelengths) and is 
not the same as used previously, 
b. Reflection. 

Radiance due to reflectance depends on the illumi- 
nating source. This source may be the sun, active, or semi- 
active sources. It is obvious that at night for. a passive 
infrared system, the radiance due to reflection is not a con- 
tributing factor. For this reason only the radiance due to self- 
emission is considered in this section. 

5 . Target Temperature 

Since the self-emittance of a target depends on the 
temperature of the target, a method for determining this tem- 
perature is needed. The temperature of an aerial target varies 
depending on the aspect of the target. The propulsion system 
has hot surfaces such as the nozzle and exhaust plumes. There 
may also be hot surfaces due to aerodynamic heating and/or 
solar radiation. 
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Air breathing engines normally have exhaust plumes 
ranging from 600 to 1000°K. Rockets typically have much hotter 
plumes. The flame temperatures for liquid propellants range 
from 2500 to 7500°K. Solid propellants flame temperatures range 
from 1700 to 3500°K. The plume temperature can be estimated 
from the relation. 



T 

flame 



T plume 




1 + Me2 



where T^ = stagnation temperature 
T^ = static temperature 

M = Mach number at the nozzle exit 
e 



a. Example I 

For a flame temperature of 2700°K and an exit Mach 
number of 3.0 the plume temperature can be found, 



T plume 



T 

flame 



1 + 



Y-l 



M 



1270. 6°K 



For many missile encounters the exhaust plume may 
be shielded from the infrared sensor. For a head-on encounter 
the temperature of interest is the skin temperature of the target. 
This temperature is due to aerodynamic heating and is a function 
of the target speed and the target material. One approach to 
finding this temperature is through the use of the recovery 
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factor, which requires some knowledge of the material of the 
target. The recovery factor, r, of a material is defined as 
follows : 



T - T 

surface ambient 

T . - T 

stagnation ambient 



The skin temperature of the target then becomes, 



T =T + r (T -T) 

surface ambient stag amb 



The stagnation temperature is found from the relation, 



T 

stag 




The Mach number, M, is that of the target and the specific heat 
ratio, y, is for air. 

b. Example II 

A target flying at M = 2.5 where the ambient temper- 
ature is 300°K has a recovery factor of 0.75. 



T - 
stag 



T surface 



575°K 



= 581°K 



This is the temperature used, along with the emissivity of the 
target, to find the radiant emittance of the target from 
equation (4) . 
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6 . Simple IR System 



Target 
& t / e, T 



IR 






Detector 



4 



filter 



chopper 



l 



preamp 



Lens 



< R > 



Figure (4-7) . Simple IR system. 



Figure (4-7) , above, shows a simple IR system and a 
target at a range R. If the system is sensitive to radiation 
in the 3 to 5 micron region the radiant emittance becomes, 



a 2 = 5 



w = 



V 3 



f e (W.) dA 
\ A BB 



5 , 

f [2vhc 1 

J La 5 exp (he/ kT) -1 J 



dA 



The above integral is best evaluated on the computer. If the 
ambient temperature and target speed is known, the skin temper- 
ature of the target can be determined. 
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The radiance from the target then becomes. 



N = 



W 

7T 



If Figure (4-7), the solid angle of the aperture as seen from 

the target is, ^ _ _a 

R 2 

The power seen at the detector surface is then, 

IR Power = N ft 

The above formula assumes no attenuation by the atmosphere or 
the IR system optics. This attenuation is significant in actual 
IR systems. These factors are normally accounted for through 
the use of atmospheric and IR optics transmission coefficients. 



T^ = transmission of the atmosphere 
T ir = transmission of the IR optics 
The total power at the detector then becomes , 



P 



T 

a 




A 

a 



W 



(5) 



7 . Detectors 

Detectors are devices which are radiation transducers. 
It's purpose is to change the incoming radiant power to an 
electrical signal, which can then be amplified. Detectors can 
be divided into two main categories. (1) Thermal detectors - 
The responsive element of a thermal detector is sensitive to 



114 



temperature changes brought about by the incident radiation. 

(2) Photodectectors - Responsive elements of photodetctors are 
sensitive to the number of incident photons. 

Detectors also are made up of windows, apertures and 
Dewar flasks. The window restricts the bandwidth to which the 
detector is sensitive. The aperture may limit the field of vie- 
in order to limit photon noise. The Dewar flask cools the 
detector which improves the detectivity. 

Detectivity of a detector is defined as, 

_ signal/noise _ S/N .... 

D = — — — ~ — \ b ) 

input power P 

The specific detectivity is 

D* = D [ Af A d ] 1/2 (7) 

Af = Bandwidth 

A. = Sensitive area of detector 
d 

For a tracking system the bandwidth is that of the preamplifier 
in Figure (4-7) . The input to the preamplifier is proportional 
to the incoming IR energy, which has been modulated to give tar- 
get resolution from the background and provide line-of-sight 
information. 

A simple chopper is shown in Figure (4-8) . It consists 
of an opaque material which has a wedge cut-out of angle a. The 
rotation causes the input from a point source to be modulated. 
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w 



Figure (4-8) . Simple Chopper. 



while that of the background is not. The input signal to the 
preamplifier would look like Figure (4-9) . The frequency content 
of the signal in Figure (4-9) can be found from a Fourier Analysis. 
If the pulses are assumed to be sinusoidal of period T, the 
optimum bandwidth is, 




T 



W 



tt aw 



s 




Background 



T 



T 



w 

s 



Figure (4-9) . Preamp input. 
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The specific detectivity is characteristic of the 



detector used. Reference [9] is an excellent source of infor- 
mation on operational detectors. 

8 . Idealized Range 
Equation (5) is, 



P 




W 



From equations (6) and (7) this becomes, 



D* 

[Af A d ] 



1/2 



S/N tt R" 



T T tt , A. A W 
a IR t a 



The idealized range, where the signal to noise ratio is unity, is 

D* T T A. A W 
„ 2 a IR t a 



o •“ “ " : 1/2 

[Af A d ] it 



To simplify this equation the radiant intensity is given by, 



J = 



A W 

TT 



Since the parameter of interest in missile design is the aperture 

ttD _ 



diameters A , it is replaced with 



so that. 



R 0 = 



D = 
a 



D* T T td tt D J 
a IR a 



-.1/2 



4 v/Af A d 

R 0 2 

D* T T td it J 
a IR 



-1 1/2 
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a. Example (Idealized Range) 

From example II a target flying at M = 2.5 has a skin 

temperature of 581°K. If this target has a presented area of 
2 

1 m , the detector size needed to detect the target at a range 
of 10,000 m can be determined. From Wien's displacement law the 
maximum radiation occurs at. 



X 

max 



4.99 \i 



If the system is designed to accept radiation from 3 to 5 microns, 
and the emissivity of the target is 0.7. 

^5 2 

TT _ _ 2-rThc 1 ^ 

£ J .5 exp (hc/XkT) -1 dA 

3 A 

If this equation is integrated on the computer, the radiant 
emittance becomes, 

W = 1410 watts/m 

The radiant intensity becomes, 

A W 

J = — — =448.82 watts 

TT 

Typical values of the parameters in the idealized range equation 
are (9) 



T =0.75 
a 

t ir - °- 95 

Af = 1000 cps 

* , 2 
A, = 1 cm 
d 

D* = 1 x 10 10 (cps) 1//2 



cm/watt 



Substituting these values into the idealized range equation gives. 



D = 0.0355 m = 1.4 in. 
a 

As can be seen the size of the IR seeker is relatively small 
compared to other seekers. The above analysis is for the 
"idealized" range. Attenuation of IR radiation can be quite 
high thereby increasing the seeker size required. 

F. WARHEAD SIZING 

The conditional kill probability of a missile is the prob- 
ability that the target is destroyed given that the warhead is 
delivered to a point in space and the fuze detonates the war- 
head at a miss distance r. 




Figure (4-10) . Encounter geometry. 
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The fragment distribution, D(<J>), is the number of fragments per 
unit solid angle, fi. The total number of fragments within the 
cone is, 

N = D (<J>) fl 



The fragment density, p, is the number of fragments per unit of 
area normal to the path. 

n D($)n 

P A A 



however, 




and 



P = 



_ D(cft) _ 



D U) sin (f> 

2 



m 



( 8 ) 



1 . Target Vulnerability 

The vulnerability of aircraft or missile components is 
normally determined experimentally. Fragments of a specified 
size are fired at the component, and the damage is assessed to 
determine if the fragment would cause a kill. If a large num- 
ber of fragments are fired, the ratio of killing fragments to 
hits can be determined. This ratio is defined as the probability 
that given a hit a kill will result, p !k _ Assuroing the 

*V H H 

distribution of hits is uniform over the target. 



K/H 



( 9 ) 
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Ap = presented area of target 



A v = vulnerable area of target 

As would be expected, depends greatly upon the 

encounter geometry. It will be dependent upon the aspect of 
the aircraft, and also depends upon the type of kill specified. 
If the target is assumed to be a spherical target, the prob- 
ability of kill given a hit can be assumed constant. 

2 . Conditional Kill Probability 

The fragment density is given by equation (8) . From 
this density expression, the average number of hits, a, on the 
vulnerable area of the target is, 

2 

D (4>) sin <J> _ 

a = — — — * — A 



It is customary to assume that the distribution of hits on the 
presented area follows a Poisson distribution. The conditional 
kill probability then becomes, 



P 



D 



P 



D 



1 

1 



-a 

e 



exp 



-D (<j>) sin 4) A v 

2 



m 



(10) 



The above expression depends upon the fragment distribution, 

D(<}>) . If the warhead casing is scored such that it produces N 
fragments of uniform size and mass, m, the problem is simplified 
by formulating an alternate expression for a. 
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£ + 2r tan 3 

w m 




Figure (4-11) . Static encounter. 

The warhead-target encounter geometry is shown in 
Figure (4-11) . As shown the area of the fragment ring depends 
on the miss distance, r^. The area of the ring is given by, A^. 

A_ = 2tt r (£ + 2r tan 3) 

f m w m 

If the N fragments are distributed evenly in A^, the number of 
fragments per unit area is 

_ N _ N 

P A. 2it r (1 i + 2r tan 3) 
f m w m 



The average number of hits on a vulnerable component is then, 
a = p A 

v 
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P K/H A p 



From equation (9) , A = 

Therefore, a = p P„ A 

K/H p 

= H P K/H A o 

a 2 7 T r (£ + 2r tan 6) 

m w m 



If in Figure (4-11) , the average target width into the paper is 
given by W, the presented area is, 

A = W (£ + 2r tan 8 ) 

p w m 



The average number of hits is then. 



_ N P K/H W 

a i “ ~ 2F~r 

m 



( 11 ) 



The distance at which the target just fills the fragment ring 
is the critical miss distance, and the maximum distance for 
which a^ applies. This is found by setting 



L = £ + 2r tan 8 

w c 

From which 

L - £ 

r = 

c 2 tan 8 



If the miss distance in Figure (4-11) is such that the entire 
target is always presented to the fragment ring; i.e., r^ > 
the presented area becomes, 

A = W L 
P 
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The average number of hits then becomes, a 



II 



a 




( 12 ) 



II 2tt r Vi + 2r tan B) 
m w m 



3 . Sizing the Warhead Radius 

The parameter of interest in this chapter is the diameter 
of the warhead required to achieve a specified kill probability. 
The development thus far is the conditional probability of kill, 
Pp. It has been assumed that the guidance system delivers the 
warhead to the point of interest, and the fuze detonates the 
warhead at this point. Since only the conditional probability 
of kill is determined, the purpose of this section will be to 
maximize or to find the warhead diameter which sets P_ = 1. 



from an alaysis of the threat, the siz e and impact velocity of 
the fragments necessary to kill the target can be determined. 
The initial velocity required to obtain the impact velocity is 
a function of the explosive used and the charge to mass ratio, 
C/M. 



D 



D 



From the threat to be encountered, the target presented 




C _ Mass of explosive/unit length 

M Mass of warhead casing/unit length 



The initial velocity is given by Gurney's equation 




( 13 ) 
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Gurney's constant, -J 2E, for various explosives are given below [11] . 



Explosive Density, kg/m 2 



TNT 


1590 


RDX 


1650 


HMX 


1840 


PETN 


1730 


Tetryl 


1620 


Composition B 


1680 


Octol 


1800 


From equation (13) , the 


charge 


the specified initial velocity 


■ c _ V i 2 / 2E 




M 1-V i 2 /2(2E) 



v/2E , m/sec 

2316.5 

2834.6 

3118.1 

2834.6 
2500.0 

2682.2 

2895.6 

to mass ratio necessary to attain 
can be determined. 




Figure (4-12) . Warhead. 
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From Figure (4-12), C/M can be expressed as. 



where 



C 

M 



2 

ur p 
e e 



, _ 2 2 , 

<* r w - *r e )p c 



P e = explosive density 
P c = casing density 



C 

M 




(r +t) 2 - 
e 



2 

r 

e 




C 

M 



2 

r 

e 



2r t + t 2 
e 




P 

c 



If t < < r 

e 



C 

M 




2r t p p 2t 
e c c 



C. 

M 





t) 



(14) 



From equation (14) the casing thickness in terms of the war- 
head radius can be determined. 



t = 



w 



Pc + 
M p 0 



1 
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With the casing thickness fixed as a function of warhead radius 

the number of fragments, N, can be determined. In order to 

achieve the desired velocities of fragments, it is essential 

to have enough length for the diameter. An acceptable length 

to diameter ratio for a cylindrical warhead is from two to three. 

Most air-to-air missiles have a length to diameter ratio of 2.5. 

For this analysis a value of 2.5 is used. 

The warhead casing volume is then given by , V , 

c 

V - 2irr 1 1 

c w w 

V = lOirtr 2 

c w 

The mass of the casing is, m c , 

m = lOirtr 2 p 
c w c 



The total number of fragments, N, is obtained by dividing the 
case mass by the individual fragment mass. 

m lOirtr 2 p 
... c w c 



From the equation (11) or (12) the average number of hits can 
be determined as a function of warhead radius 



a l 



5tr P P T . ,„W 
w c K/H 



mr 



m 



5tr p Pr, /tl WL 
_ w c K/H 

11 mr U +2r tang) 
ro w m 
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The conditional probability of kill can be determined now for a 
given warhead radius using equation (10) . 



a. Design Example (Effect of radius on P^) 

If the above equations are programmed for the con- 
ditions listed below, a plot of the conditional probability of 
kill versus warhead radius can be obtained. 

Threat: RAM-K 



L = 19.51 m 
w = 6.10 m 



P K/H 



.10 



Fragments : 

m = 105 grains = 0.0068 kg 
= 2133.6 m/sec 
P c = 7000 kg/m 
6 = 20 degrees 



Explosive: Composition B 

P e = 1680 kg/m 

Miss distance: r = 50 ft = 15.24 m 

m 



Figure (4-13) is a plot of the output. From this 

figure it can be seen that a warhead radius of r = .06 m is 

w 

required to achieve a conditional kill probability, P D = 1.0. 
Therefore the missile diameter required for warhead considera- 
tions is, d = 4.72 inches. 

The warhead radius also varies with the required 
initial velocity. From the equations for charge to mass ratio 
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(meters ) 



and casing thickness the radius required to achieve a specified 
kill probability can be determined. 



Since, 



and 



t = 



r 

w 




c = V i 2 / 2E 

M 1 - V i 2 /2(2E) 



the average number of hits becomes, 



a l = 



5 P P v /ti w r 3 
' c K/H w 



V i 2 /2E 



-£■ + 1 



,1 - V ± /2(2E )/ M e 



mr 



m 



The above equation assumes r < r 

me 



Letting 



b = 



5 c P K/H W 



'V. 2 /2E 



,1 - V ± 72(2E), 



— + 1 



mr 



m 



the probability of kill becomes, 

2 



-br 



P D = l-e 



W 



If the conditional probability of kill is selected as P 



equation (15) can be solved for the warhead radius. 



(15) 

= 0.999, 



130 



r 



w 




(16) 



The initial velocity, V. , required to achieve a target kill is 
a function of the miss distance, fragment size, expected en- 
counter altitude and target characteristics. For this simple 
analysis the effect of initial velocity on warhead radius will 
be studied. The effect of varying the initial velocity in 
equation (16) is plotted in Figure (4-14) . The results were 
determined for various explosives to show their effect on the 
warhead radius. If an initial velocity of 7000 ft/sec is chosen 
with Composition B as the explosive the required warhead radius 
is 2.6 inches. 

From this analysis the warhead weight can be found. The 
casing thickness is given by. 



t = .115 r (from required C/M) 



w 



l = 5 r 



w 



w 



Therefore, the explosive weight is 



W = tt (r 
e 



w 




W = 5.9530 kg 
e 



The casing weight becomes 



W = (irr 
c 



2 




w 



W =6.8651 kg 
c 
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Figure (4-14). Effect of initial velocity on warhead radius. 



The total warhead weight becomes, 



W T7 „ = 12.8181 kg = 28.26 Jib 
WH 3 

The required diameter for the radar antenna was 10 
and exceeds that required for the warhead. If the warhead 
made hollow and kept at the same weight, equation (14) can 
reformulated to give 



C 

M 



l"(r - 


, 2 


21 


t) - 


r. p 


L w 


i J e 


T2r t 


- t 2 " 


P 


L w 




c 



The hollow portion of the warhead is then found from 




) 



ilw 



Figure (4-15) is the resulting warhead. 




Figure (4-15) . Hollow warhead. 



inches 

is 

be 
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V. BASELINE DEFINITION 



A. CONTROL CONCEPTS 

The lifting and control surfaces of a missile may be of 
monowing, triwing or cruciform configuration. Figure (5-1) shows 
these three arrangements . 




Monowing Triwing Cruciform 

Figure (5-1). Control configurations. 

The monowing arrangement is typical of most cruise missiles, 
which require long range and low drag. For this type of arrange- 
ment the missile must bank to orient the lift vector for a maneuver. 
Because of this, the monowing missile is not as rapid in maneu- 
vering as the cruciform configuration, which can produce lift 
in any direction instantaneously. The cruciform control also 
has identical pitch and yaw characteristics which results in 
a simpler control system. The triwing configuration is used 
very seldom for conventional missiles. It can be shown that 
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the triwing missile requires larger wing size; therefore, there 
is little drag savings even though there is one less wing 
than with the cruciform missile [10], From this discussion it 
can be seen why the cruciform configuration is the most com- 
monly used for tactical missiles of short or medium range. 

The positions of the lifting surfaces on the missile body 
depends on the method of control used for the missile. There 
are three conventional methods of control for tactical missiles. 
These are; 1) Wing control; 2) Tail control; 3) Canard 
control . 

1 . Wing Control 

Wing control missiles normally have large movable wings 
located slightly behind the missile center of gravity. Because 
of the small moment arm the wing surface must be relatively 
large to provide control effectiveness. As would be expected, 
larger activators are required for moving the wings. A positive 
deflection of the wings causes a positive normal force; there- 
fore, the missile reacts almost instantaneously; thus making 
wing control the fastest reacting method of control. Because of 
the smaller moment arm, the resulting smaller pitching moment 
and instantaneous lift result in smaller angles of attack. This 
feature makes wing control missiles attractive for applications 
where theincidence angle must be kept small. Air breathing 
applications often use wing control because of inlet perfor- 
mance degradation at higher angles of attack. This type of 
control is also good for fixed seekers. 
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2. Canard Control 



Canard control is normally the method of control where 
the movable surfaces are placed as far forward as possible. 

Because of the resulting larger moment arm, smaller surfaces are 
required to provide control effectiveness. Lift on the missile 
is still developed primarily by the aft (wing) lifting surface. 
Response is slower than wing control because of the need to 
pitch the missile to an angle of attack before lift is developed. 
Higher angles of attack are needed to generate the required lift. 

One advantage of this type of control is convenience of packaging. 
Since the controls and avionics are forward of the propulsion 
system, the need for connectors is elminated. For stability 
reasons the wing of a canard control missile must be located 
farther aft than a conventional wing-tail configuration. The 
zero lift drag is normally lower than wing control missiles 
due to smaller surfaces. 

3 . Tail Control 

Like the case of canard control, a tail control missile 
has the movable surface as far from the center of gravity as 
possible. This also results in a larger moment arm and there- 
fore smaller surface required. Control deflection is the oppo- 
site of that for canard or wing control since a negative control 
deflection results in a pitching moment that pitches the nose 
up and therefore a positive lift on the main lifting surface (wing) . 
Tail control is normally the slowest method of control. One 
advantage is that the flight controls are at the end of the 
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missile requiring the propulsion system to be located forward 
of other types of control. This results in less center of gravity 
movement as the propellant grain burns. This becomes a definite 
advantage for longer range missiles. 

B. GROSS WEIGHT AND CENTER OF GRAVITY 

Since the component weights and their precise locations can- 
not be determined at this point in the design process, some 
method of estimating the total weight and center of gravity is 
needed. There are two approaches commonly used to find the gross 
weight of missiles in the conceptual design phase. One is through 
the use of the historical data discussed in Chapter 2. Since 
the warhead weight is known, an estimate of the gross weight 
is now. 



Another method, which is used extensively in the design of 
aircraft, is the use of regression formulas to find gross weight 
or component weights in terms of parameters that are known early 
in the design. Reference (7) has derived such a formula for the 
gross weight of a missile. It is, 




w G = K g (L) 



2.13 



(D) 



1.14 



K = Constant to be determined 



( 1 ) 



L = Total missile length (inches) 



D = Missile diameter (inches) 
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The constant K in equation (1) is derived from a baseline 
G 

(generic) missile. In this case 



K G = 



(L) 



( w G ) 



2.13 



baseline 



(D) 



baseline 



1.14 

baseline 



The accuracy of a regression formula such as equation (1) depends 
upon how close the synthesized missile is to the baseline for 
which K was determined. If the parameters used vary more than 
20 percent, the accuracy of the regression equation decreases 
rapidly; therefore, if the parameters length and diameter vary 
significantly from those of the baseline, equation (1) may give 
an erroneous estimate of the gross weight. As an aid in deter- 
mining the gross weight, the following values of K were derived 
from data given in reference [7] . 



Missile 



Figure 



K, 



SRAAM 

MRAAM 

LRAAM 

SAM 



(5-2) 

(5-3) 

(5-4) 

(5-5) 



.00128 

.00118 

.00093 

.00108 



Figures (5-2) , (5-3) , (5-4) , and (5-5) show the generic missiles 

from which these values were derived. 

The center of gravity of the baseline missile for this section 
is taken to be 60 percent of the total length. At this point 
sufficient information has been developed to define the baseline 
missile from which design iterations will be made. The lifting 
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Figure (5-3). Generic MRAAM configuration [7]. 
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Figure (5-5). Generic SAM configuration [7]. 



surface planform for the baseline is taken as a delta planform. 

As will be shown later, this wing planform will be very close to 
an optimum wing. 

C. DESIGN EXAMPLE (BASELINE DEFINITION) 

The threat for which this missile is designed is a highly 
maneuverable fighter. For this threat a canard control, cruci- 
form configuration is chosen. The rationale is that this configur- 
ation will provide the fast response necessary at minimum drag. 

The diameter was fixed at 10 inches due to antenna considerations 
in Chapter 3. From historical data in Chapter 2, 

L = (^-) D = (15.89) 10 = 158.9 inches 
L 

L^ = (— ^) D = 22.3 inches 



The MRAAM of Figure (5-3) is used as a generic missile for 
selecting K_. Inserting length and diameter into equation (1) , 
the first estimate of gross weight becomes, 

W„ = .00118 (158. 9) 2,13 (10) 1 * 14 = 794.83 lb 



The total lifting surface required is then, from historical 
average, 



S ^AVG (Wq) ( 88.09 



) (794.83) 



S = 9.02 ft' 
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Since the canard to wing area ratio is known from historical 
data , 



S /S„ =0.20 
c W 

S r7 = 7.52 ft 2 
W 

S = S - S T7 = 1.50 ft 2 
c W 



Lifting surfaces (delta planform) . From historical data. 



AR W - !- 61 

AR =3.74 
c 




AR 




s 



Figure (5-6). Lifting surface. 



Wing: (b) w = \f AR W S w = 3.48 ft 

(C ) = 2 S /b = 4.32 ft 

W 

Canard : 

(b) =/AR S = 2.37 ft. 
c c c 

(C ) = 1.27 ft. 

c 
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The baseline missile is now defined. The canards are placed 
as far forward as possible. The wings are placed as far aft 
as possible to ensure the center of pressure is behind the 
center of gravity. The exact location of the wings will be 
modified in the next chapter. Figure (5-7) is a drawing of 
the baseline missile. 
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VI. LINEAR AERODYNAMICS 



The total drag on the missile is used to size the propul- 
sion needed for the cruise segment of flight. The wing size 
depends upon the maximum lift required by the missile. The 
wing and tail are normally placed to provide minimum drag during 
cruise or a certain stability margin at launch. To make any 
of the above calculations, values of the aerodynamic coefficients 
are needed. This chapter presents the background aerodynamic 
theory necessary for these initial calculations. The theory 
used is linear aerodynamic theory and slender body theory, from 
which simple calculations can be made. Where linear theory did 
not apply, an attempt was made to find existing empirical expre- 
sions, which yield results accurate enough for initial calcula- 
tions. The full nonlinear theory will be presented in Chapter 9. 
The reference area for all coefficients in this report is the 
missile maximum cross sectional area. The reference length is 
the maximum missile diameter. 



A. MISSILE DRAG 

The total drag of a missile consists of zero lift drag, C , 

D 0 



and induced drag, C 



D. 



C = C_ + c 
D D 0 D i 



( 1 ) 



The zero lift drag can be found from a component build up method 
in which the contributions due to the nose, body and lifting 
surfaces are added together to obtain the total zero lift drag. 
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Care must be taken to reference the appropriate areas when 
using this method. The total is then multiplied by 1.25 to 
account for interference effects and variations in skin roughness 
[ 12 ] . 



'D, 



= 1.25 



‘S' 



N 



( V 



B 



‘S’ 



w 



V 



( 2 ) 



The method used to find the components drag depends upon the 

speed regime in which the missile is operating. Since this 

report is concerned primarily with supersonic tactical missiles, 

supersonic zero lift drag will be discussed here. The component 

supersonic zero lift drag can be divided into skin friction, 

C , and wave drag, C . 

D f D W 



(C ) = <C > + <c ) 

0 w f w w w 



1. Supersonic Skin Friction 

The flow over a body traveling at supersonic speeds is 
likely to be turbulent; so the incompressible skin friction co- 
efficient is given by, 



•f . 
i 



(log 



.455 

10 V 



2.58 



(3) 



The Reynolds number in equation (3) is based upon the cruise 
altitude and speed and upon the characteristic length for the 
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component being determined. The Reynolds number is given by, 



R e = 



p V M X 



where x = , the length of the nose 

L , the length of the body without the nose 

(c) , the mean aerodynamic chord 



The compressibility correction to equation (3) is, 



-0.58 



C f = C f (1+0.15 M M ) 

l 



(4) 



From equation (4) the skin friction drag coefficient for each 
component can be found when referenced to the appropriate area, 

(C ) = (C ) wetted } N 

f N N ref 



^wetted^ 



< C D > - (C £> S 7" 

f E B ref 



B 



The lifting surface skin friction drag is determined in a similar 
manner. Care must be taken to include all surfaces in the 
wetted area calculation for the lifting surfaces. 

2 . Supersonic Wave Drag 

The supersonic wave drag consists of components contri- 
buted by the nose and lifting surface. Nose wave drag depends 
on the shape of the nose, and the most common nose shapes are 
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conical, ogival and hemispherical. Reference [10] lists 
empirical formulas for finding the form (wave) drag of various 
nose shapes at zero angle of attack. 

(1) Conical; C D = (0.083 + 0 . 09 7/M 2 ) (a/10 ) 1 * 69 

u w 

a = tan -rp— = Nose semi-vertex angle (Degs) 

2L N 

The center of pressure for a conical nose is at the centroid of 
the nose planform or two thirds the length of the nose. 

( 2[196(L /D) 2 - 16] 

(2) Ogival; C n = P< 1 ~ 

W ( 28 (M + 18) (L N /D) 

P = (C_ ) for conical nose 
D W 

The center of pressure for an ogive noise is, 



fp = 1 50 (M + 18) + 7M 2 P ( 5M - 18) 
L N 2 40 (M + 18) + 7M 2 P (4M - 3) 



The semi-vertix angle for an ogive is twice the equivalent cone 
angle. 

(3) Hemispherical; The drag on a hemispherical nose is extremely 
high compared to other nose shapes, and is difficult to 
estimate. An initial estimate of the wave drag can be 
found from Figure (6-1) . 
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Figure (6-1) . Hemispherical nose wave drag. 



The wave drag due to the lifting surfaces can be found 
using the methods of reference [12]. For a double wedge air- 
foil with sharp leading edges as shown in Figure (6-2) , the 
wave drag is given by the following formulas: 



(1) Supersonic leading edge. 



W 



= | (§) 
6 c 



W 



ref 



(2) Subsonic leading edge. 



2 S 

C D = B COt A le (§> g-S- 
W ref 




<r 

<r 




■> 




Figure (6-2). Double wedge wing [12]. 



where „ _ C//x t 

= 1 - x t /c 




A^ E = leading edge sweep 
S w = planform area 
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3. Base Pressure Drag 



The drag contribution of the blunt base for non-boattailed 
bodies is given in reference [4] as. 




This term assumes no jet thrust from the base of the missile, 
or that the missile is operating in the power off condition. 
This term is not included in equation (2) , which is not a bad 
assumption for powered flight where the nozzle exit area is 
approximately equal to the base area of the missile. For the 
case where the nozzle exit area is much less than the base 
area as in Figure (3) , the base pressure contribution should be 
included. 





Figure (6-3). Base pressure areas [10]. 
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In this case the base drag is, 



C 



AB 



C 




= shaded area of Figure (6-3) 



4 . Induced Drag 

The induced drag on a missile is the drag due to lift. 
This drag is caused by the component of the lift vector in the 
drag direction. For supersonic flow the induced drag is given 
by. 



The lift curve slope, C , will be developed later in this chapter. 



B. DESIGN EXAMPLE (ZERO-LIFT DRAG CALCULATION) 

The thickness to chord ratio of the wing and tail have not 
yet been determined; however, it is desirable to construct the 
lifting surface as thin as structurally possible to minimize the 
wave drag. Since structures have not been covered the minimum 
thickness to chord ratio is estimated at 3 percent. The flight 
and geometric conditions determined thus far are; 



C 




where 



W 



c 



L q S j. 
^ ref 
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h = 10,000 ft 


L rT = 1.8583 
N 


O 

• 

ro 

II 

s 


L cb = 11.3833 


S . = 0.5454 
ref 


(c) w = 2.8816 


(t/c) w = ( t/c ) c s 0 . 0 3 


(5) = 0.8467 


S T7 = 7.52 ft 
W 

S c = 1.50 ft 





From these conditions the following table can be constructed. 





Nose 


Afterbody 


Wings 


Canards 


X 


1.8583 


11.3833 


2.8816 


0.8467 


S wetted 


3.3394 


29.8014 


15.0400 


3.0000 


R 

e 


2.984 


18.282 


4.6280 


1.360 


c fi 


0.0025 


0.0020 


0.0024 


0.0029 


c f 


0.. 0015 


0.0012 


0.0015 


0.0017 




0.0095 


0.0652 


0.0827 


0.0187 


C ^W 


0.1545 


— 


0.0400 


0.0050 


O 

Q 

U 


0.1640 


0.0652 


0.1227 


0.0237 


The zero 


lift drag 


for the wings 


and canards 


in the above table 


take into 


account 


that there are 


two sets of 


wings (4 panels). 


The total 


missile 


zero lift drag 


then becomes 


t 



= 1.25 (C ) + (C ) + (C ) 

u o u o u o 

U N U AB U W 



+ (C_ ) 



C D = 0.4695 



155 



C. MISSILE LIFT CURVE SLOPE 




Figure (6-4) . Wing-body- tail lift [13] . 



For the purposes of this section lift on a wing-body-tail 
combination can be taken as the sum of the components in Figure 
(6-4). These consist of. 



L. t = lift on the nose 
N 

L. = lift on the wing in the presence of the body 
W (B) 

Lg = additional lift on the body due to the presence 
(W) of the wing 

L = lift on the tail in the presence of the body 

X (B) 

L = additional lift on the body due to the presence 

B (T) of the tail 



The lift of only the wing-body combination can be defined as, L c , 
where , 
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The lift of the wing alone, L^, is that obtained from thin 
airfoil theory or experiment considering the exposed area only. 
The constant K c in equation (5) is defined as 



K _ = K„ 
C B 



+ K, 



(W) 



w 



(B) 



From this equation, K is the ratio of the additional 

B (W) 

body lift in the presence of the wing to the wing alone lift for 
zero control deflection, 6=0. 



B 



(C. ) 
La 



K 



(W) _ 



B 



(W) 



B (W) L W 



^ C La^ 



W 



W 



K. 



(B) _ 



(C La’ 



W 



W 



(B) 



W 






(B) 



W 



6 = 0 



6 = 0 



The interference factors, K and K_ T have been deter- 

B (W) ™(B) 

mined from slender body theory for wing-body combinations. 

Figure (6-5) is a plot of these values. The wing alone lift 

curve slope, (C T ) is determined from thin airfoil theory or 

La W 

experiment. The lift curve slope of delta wings with super- 
sonic leading edges is given by. 



C Lct = 6 wher ® 6 = 



"m 2 - 1 



k < e 



For subsonic leading edges this becomes. 



'La 



2irtan£/tanu 

E6 



k > e 
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Figure (6-5). Interference factors [12], 
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Where E is the elliptic integral of the second kind for 



'J 1 - ( tane/tany) , 

Figure (6) shows the Mach angle y and sweep angle, e. 




Figure (6-6). Wing leading edge. 

The lift of the nose is that obtained from slender body 
theory. For small angles of attack 

(C Na> “ 2/rad - 
N 

The tail-body combination lift is determined in the same 
manner as the wing. The above equations are defined for the case 
where the incidence angle is zero and the angle of attack, (a). 
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is varied. Analogous terms can be defined for the case where 
the angle of attack is zero and the incidence angle (6) is varied. 



B 



(C L6 ) 



(T) _ 



B 



^B 



(T) 



"^7 



(T) 



(B) 



(C L6> 



(B) 



(B) 



(C La ) 



a = 0 



a = 0 



The interference factors k and k were also found 

B (T) a (B) 

from slender body theory and are plotted in Figure (6.5). 

From the above defnitions the total missile lift curve slope 
can be found. If the analysis is for small angles of attack 
so that L~N. 



(N) cm = (N) N + (N) W + (N) B + (N) T + (N) 
Cm N W (B) B (W) X (B) B 



(T) 



in coefficient form, 



C N q S reg = <C N> <1 S ref + (C N> q S w + (C N> q S W 

N ”(B) B (W) 



(C N>, 



g S + 
4 t ° T 



(C N» 



(B) ) 



B 



g S 
4 t 



(T) 



If the downwash is neglected so that q T = q , and the above 
equation is differentiated with respect to angle of attack. 
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s s 

(c No > - (Sa 1 + ’Sla’ + 'Sa 1 

a CM Na N a W,_. ref a B .... ref 

(B) (W) 



+ tC Na* S . + (C NC? S 

T (B) ref B (T) ref 



From the definition of the interference factors 



(C XT J = (c„j + (c 



'N of 



CM 



Na' Na) 

N 



V., + k b,„ > 



w 



W (B) °(W) ref 



+ (C Na } (K T + K B } S“^ 
Na T 1 (B) B (T) ref 



( 6 ) 



If the tail is the control surface, a similar development for 
the control effectiveness, ^ yields x 



(C N6 } = (C Na ) T (k T + k B } s - ^ 

N6 CM Na T T (b) B (t) S ref 



D. MISSILE PITCHING MOMENT 

From Figure (6-7) the moment about the center of gravity 
can be found. The centroid of the wing is the location of the 
center of pressure of a wing alone in supersonic flow. The 
effect of the wing-body combination is to move the center of 
pressure aft and as can be seen in Figure (6-7) the center of 
pressure for the additional lift on the body due to the presence 
of the wing is aft of this location. Reference [13] has an 
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B (W) 



x, 



T (B) 



x 



B (T) 



Figure (6-7) . Forces Acting on the Missile. 



excellent discussion on how to find these centers of pressure. 

For the purpose of this chapter, which is initial sizing and 

placement of the lifting surfaces, the center of pressure for 

both of these forces, N rT and N„ is taken as the centroid 

"(B) B (W) 

of the wing planform. With this assumption, the moment about 
the center of gravity is. 



(M) = N m x m + (N T7 

CM N N W 



(B) 



+ N fi > *W + (N T_ + N B > X ' 



(W) 



(B) 



(T) 



Following the same development as for C 



N 



, X N, S - 



C Ma = (C Na> -5 + (C Na> (K W, B1 + *b > ‘TT* <r*I> 

N W (B) (W) ref 



+ (K T, > + K B, ( s 

T (B) (T) ref 



(7) 
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Also if the tail is the control surface. 



'M<5 



= (C ) (k + k ) (5jf> (g^- 

Na T T (B) B (T) D S ref 



) 



Ome must be careful in defining the moment arms in the 
above equations. If a nose up pitching moment is developed in 
Figure (6-7), the moment arm is positive. Conversely, a nega- 
tive moment arm means a nose down pitching moment is developed. 

E. WEIGHT AND CENTER OF GRAVITY VARIATIONS 

From the preceding sections of this chapter, it can be seen 
that the analysis is normally for one point in the flight pro- 
file. Since the launch condition has been defined, the tail 
sizing can be accomplished for this initial condition. The 
variation of missile weight and center of gravity is due to 
propellant burning. For a solid propellant this variation can 
be quite large. Since no information is available on the pro- 
pellant at this point in the design, the following guidelines 
will serve for initial calculations. For air-to-air missiles 
the propellant is approximately 35 percent of the launch weight. 
The center of gravity travel is 5 percent of the length. For 
surface-to-air missiles which must be boosted to flight speed, 
the propellant weight can be taken to be 48 percent of the 
launch weight with a corresponding center of gravity travel of 
approximately 8 percent of the length. These are initial approxi- 
mations taken from historical data and can be refined later in 
the design process. 
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VII. LIFTING SURFACE DESIGN 



From the analysis of wing lift and drag in Chapter 6, it can 
be seen that the performance of a wing will vary with planform. 

Up to this point the lifting surfaces have been considered delta 
planforms with zero taper ratio. This is not necessarily the 
optimum planform since it was taken from a historical average 
and does not apply to a specific missile. This chapter deals with 
the sizing, placement and planform definition of the missile . 



A . WING PLACEMENT 

The wing placement on the missile depends upon the type of 
control used. For canard control the wing (aft surface) is nor- 
mally fixed as far aft as possible for stability purposes. For 
a wing control or tail control missile, where the wing is near 
the center of gravity, the wing placement becomes more critical 
and depends upon the stability margin required. Since the drag 
during the cruise segment includes the drag due to lift, one 
method to minimize drag would be to place the wing such that 
zero lift is produced on the tail during cruise. The moment about 
the center of gravity is zero for trimmed flight, thus equation 
(7) yields, for (C N ) = 0. 



_ (C Na S ref 

*W (C T ) S TT (K„ +~ K 



' Na ' W W W (B) 



B 



( w) 
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B. MANEUVER LOAD FACTOR 



Regardless of the type of missile being designed, it will be 
required to maneuver in order to intercept its intended target. 
For an air-to-air or surface-to-air missile the level of this 
maneuver may be quite large. The maneuver load factor required 
was found in Chapter 3. The maximum maneuver the missile can 
sustain depends upon the maximum trimmed normal force the missile 
can develop. 



Figure (7-1) . Sustained maneuver of a missile. 

From Figure (7-1) the force developed by a missile in a con- 
stant acceleration turn is 

L = nW 



L = nW 




V 

w 



n = maneuver load factor 



( 1 ) 



Small angles of attack are assumed for which the lift and 
normal force are approximately equal. This approximation is 
good for angles of attack of up to 10 degrees, and above this 
value of maximum trimmed angle of attack the linear theory 
becomes inaccurate. With this assumption for now, equation (1) 
becomes , 



L 



N = 
(C 



Na 



(C. 



Na 



nW 

) a q S 



ref 



) 

required 



nW 



(a) 



nW 

max 



q s 



ref 



( 2 ) 



Equation (2) gives the lift-curve slope required to develop 
the required normal force at a trimmed angle of attack, ( a ^ max * 

The lift curve slope developed is given by equation (6) in Chapter 
6 . 



(C Na } = (C isra ) + (c m> (K w + K r } 

N CM N N N W "(B) B (W) S ref 



+ (C W (K T + k b } 

cga T t (B) b (T) ref 



(3) 



As discussed earlier the values of the interference factors 

depend upon the wing planform. Typical values of (IC 7 + K ) 

™(B) B (W) 

and (K + K ) are from 1.5 to 2.0. As a conservative 

(B) b (T) 

estimate a value of 1.5 is assumed. Equation (3) then becomes, 



s s 

(C. T ) = (C XT ) +1.5 (C. T ) ^ — +1.5 (C. T ) 

Na Na . T Na r7 S , Na m S , 

CM N W ref T ref 



(4) 
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Selecting supersonic leading edges. 



< C Na> w - < C Na> T ■ ^ 



Therefore 



(C Na^ 



CM 



+ 3 S ^ S W + S T^ 
N ref 



If the complete missile lift-curve slope is set equal to the 
required lift-curve slope, the lifting surface area required to 
maintain the maneuver load factor can be found, using equation (4) . 



3 S 



(S 



w 



V 



req 



■ t(C Na> 



req 



(C Na ) 



ref 



N 



Let 



(S W + V 



K = 



(S W + V 



req 



baseline 



If the same ratio of wing to control surface area is kept to 
minimize stability perturbations, the new wing and tail area are. 



S 



W 



S 



T 



K( V 



K(S t ) 



baseline 

baseline 



The above analysis assumes a linear variation of C. T which is 
good for small angles of attack. A more precise analysis will 
be performed in Chapter 9. The above analysis should be performed 
at the expected encounter conditions. 
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1 . Design Example (Maneuver Load Factor) 

Since a canard control was chosen for the design example 
the canard and wing position are fixed as fqr forward and aft 
as possible. The analysis is done for a conservative Mach number 
of 2.5 after the missile has slowed from its cruise velocity due 
to the maneuver. From previous results, 

n = 31 g‘s 

(S + S ) = 9.02 ft 2 

baseline 

W_ = 794.83 lbs 

If the missile is required to maneuver at one-half its powered 
range, and the propellant weight is 35 percent of the launch 
weight, the maneuver weight becomes 

W = 655.73 lbs 

The dynamic pressure for M^ = 2.5 and at 10,000 feet altitude is 




q = 6369.84 lbs ft 2 

The required lift-curve slope becomes for ( a max ) = 10 degrees 
= .1745 rad. 

(C_. ) = 33.5311 per rad. 

Not c 

req 

The required lifting surface area is? 

$S - 

(S TT + S„) = [33.5311-2J = 6.57 ft 

W C b 
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Therefore 



K = 0.7284 . 



The required wing and canard sizes to achieve the maneuver are; 

S rT = (7.52) ( .7284) = 5.48 ft 2 
w 

S c = (1.50) ( .7284) = 1.09 ft 2 



C. TAIL SIZING 

For a wing-tail combination, the primary concern in the 
sizing of the tail is the static stability of the missile. The 
missile becomes more stable as the mission proceeds, since the 
center of gravity moves forward as the propellant burns. As 
the missile becomes more stable, control of the missile becomes 
sensitive. If the tail is sized for zero static stability at 
launch, the missile control will remain more effective during 
flight. This is the best condition possible without the use of 
some form of stability augmentation system at launch. Therefore, 
at launch C„ = 0 . 



0 - (C Na> X N + ,C N a> *w rrr <K W. r . + K B, > 
N W ref (B) (W) 



(5) 



+ (C Na } T X T S 



-±- (K + K ) 

ref T (B) B (T) 



With the lifting surface area fixed due to the maneuver load 
factor, the tail can be sized to satisfy equation (5) . As can 
be seen the position depends highly on the moment arms that the 
lift forces act through. Since the missile length has not been 
fixed at this point in the design and may vary due to propulsion 
requirements, this analysis will be completed later. 



D. WING PLANFORM 



The wing planform is specified by the leading edge sweep, 
A tt ,, taper ratio. A, aspect ratio, AR and planform area S rT . The 
planform area was fixed due to maneuvering requirements in a 
previous section. This section is concerned with defining the 
remaining planform parameters. Figure (7-2) is the wing plan- 
form and the equations used to define these parameters. 




b 2 /s 

°t /c r 
2 s 

b (1+A) 



1 c 

3 r ' 1+A 
2 tan ^LE 



Figure (7-2) . Planform geometry. 



The value of the lift-curve slope used previously was derived 
from linear theory and is applicable only to simple planforms. 

For the purpose of this section, which will include low aspect 
ratio wings with supersonic and subsonic leading edges. Figures 
(7-3) and (7-4) are used. These figures have been corrected for 
3-D effects. 

1 . Effect of Taper Ratio and Leading Edge Sweep 

The lift and drag characteristics of the wing are the 
primary parameters of interest. The objective in wing design is 
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Figure (7-3) . Theoretical wing lift curve slope [12] . 



theory 
(per rad) 




171 



TAN 



TAN 





Figure (7-4) . Lift curve slope continued [12] . 
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is to obtain maximum lift with minimum drag. It will be shown 
that these objectives are usually conflicting; therefore, some 
compromise, or optimum, planform must be found. From Chapter 6 
the equations for the drag of a wing in supersonic flow are found. 
The skin friction drag depends upon the mean geometric chord. 
Equation (3) in Chapter 6 indicates a larger mean geometric 
chord would result in less skin friction drag for a fixed area. 

The mean geometric chord can be expressed in terms of taper 
ratio as follows, 



This equation leads to a zero taper ratio to maximize c and 
reduce skin friction if the surface area and aspect ratio were 
fixed. The wave drag is constant for supersonic leading edges, 
and decreases when the leading edge goes subsonic or the leading 
edge is behind the mach line. The lift capabilities also de- 
crease as the leading edge goes subsonic. 

As stated earlier this leads to conflicting performance since 
the objective is to minimize drag while maximizing lift. At 
this point an example best illustrates the results of varying 
taper ratio and leading edge sweep. The lift - curve slope is 
derived using the methods of reference (1). Figures (7-3) and 
(7-4) are from reference [12] and are used to find the supersonic 
linear lift corrected for 3-D flow effects. The drag methods 
of Chapter 5 are used to determine the drag characteristics. 



- _ 2 2 s 

C 3 b (1+A) 





2 . Example 



2 

For a fixed wing area of 4 ft and an aspect ratio of 
2 the lift curve slope and drag were determined for a flight 
Mach number of M= 2. This Mach number corresponds to a sonic 
leading edge sweep of 60°. The wing leading edge sweep was 
varied from zero to 75° for taper ratios of 0, 1/2 and 1. The 
results are plotted in Figures (7-5) and (7-6) . 

From this example certain generalizations can be made. 
From Figure (7-5) it can be seen that the drag is relatively 
insensitive to taper ratio. There is a reduction in drag for 
increased leading edge sweep. For this example there is approxi- 
mately a 6 percent drag reduction for every five degree increase 
beyond the sonic value. Figure (7-6) indicates that the effect 
of decreasing the taper ratio is to delay the drop in the lift 
curve slope of the wing. From this example a general guideline 
would be to fix the leading edge sweep 5 degrees beyond the sonic 
value and the taper ratio at zero. This would provide a 6 per- 
cent reduction of wave drag while maintaining the maximum lifting 
capabilities of the wing. 

3 . Effect of Varying Aspect Ratio 

The result of increasing aspect ratio is an increase in 
the lift-curve slope of the wing [12] . The aspect ratio is 
given by, 

AR = b 2 /s 

For a zero taper ratio wing this becomes? 

AR = 2b/c r 



( 6 ) 
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Figure (7-5). Effect of X on drag. 
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From equation (6) it can be seen that increasing aspect ratio 
results in a corresponding increase in wing span for a given 
root chord. Missiles are normally span limited, due to launcher 
constraints; therefore, there is a maximum span which can be 
accomodated. The drag of the wing again is in conflict with 
the lift since increasing aspect ratio decreases the mean geo- 
metric chord for a given span and therefore increases skin fric- 
tion drag. A compromise wing AR must be found which considers 
both the lift and the drag characteristics. For a discussion 
of optimum aspect ratio, the following functions are defined 
for convenience: 



where 



F = F^ + F£ ~ Lift-drag function 

F.. = jr ; — y ^Normalized drag function 

0 max 



(7) 



1/C 



F„ = 



Na 



2 ' 1/(C Na> . 

mm 



= Normalized lift function 



From equation (7) if F is plotted over the allowable 

range of aspect ratio, a minimum value of the lift-drag function 

fixes the desired aspect ratio. This aspect ratio corresponds 

to minimizing the drag function, F^ , while maximizing the lift 

function, F«. For convenience the abscissa is plotted as c/c , 

2 (ma 

as shown in Figure (7-7) . 

From Figure (7-7) the optimum mean geometric chord cor- 
responds to point A. 
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Figure (7-7) . Lift-drag function. 



c = 




c 

max 



The aspect ratio is then given by 



c 



r 



3 

2 



c 



b = 2s/c 

r 

AR = b 2 /s 



As mentioned earlier a missile is normally span limited. The 
plot of the lift-drag function F is normally fairly flat on the 
left, or for c / c ( max ) approaching zero. For this reason point 
B may be chosen as the optimum since F varies very little up 
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to this region. Point B corresponds to increased chord and 
decreased span. 

4 . Design Example (Wing Planform) 

From the previous analysis the following parameters were 

defined. 



Mm - 3-0 

S„ - 5.48 ft 2 

w 

If the leading edge sweep is fixed 5 degrees behind the mach 
line . 

y = sin" 1 ^ = 19.5° 

A LE = 95 - U = 75.5° 

X = 0 



The planform table becomes 



AR 


b 


c 

r 


2.0 


3.31 


3.31 


1.5 


2.87 


3.82 


O 

• 

rH 


2.34 


4.68 


.75 


2.03 


5.40 


.50 


1.66 


6.62 


ID 

CN 

• 


1.17 


9.36 



c 


(B) 


k b (W) 


2.21 


1.17 


.26 


2.55 


1.19 


.30 


3.12 


1.22 


.35 


3.60 


1.24 


.40 


4.41 


1.29 


.48 


6.24 


1.36 


.62 



The wave drag is constant and is given by, 

2 S 

C = 4 cot A (|) g-2- 

D W LE C S ref 



C = .0094 
U W 
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AR tan A 



6 _ 

tan A 



AR 





LE 



le 



(C„ ) (K.. +K n 

Na W "(B) B (W) 



2.0 


.0303 


.0397 


7.73 


.73 




2.1137 


1.5 


.0296 


.0390 


5.80 


.73 




1.9845 


1.0 


.0288 


.0382 


3.87 


.53 




1.8271 


.75 


.0282 


.0376 


2.90 


.73 




1.7813 


.50 


.0274 


. 0368 


1.93 


.73 




1.4809 


.25 


.0260 


.0354 


. 97 


.73 




.7937 




AR 


c/c 
' max 


F 1 


F 2 


F 






2.0 


.35 


1.0000 


.3755 


1.3755 






1.5 


.41 


.9824 


.3999 


1.3823 






1.0 


.50 


.9622 


.4344 


1.3966 






.75 


.58 


.9471 


.4456 


1.3927 






.50 


.70 


.9270 


.5360 


1.4630 






.25 


1.00 


.8917 


1.0000 


1.8917 






From Figure (7-8) 


it can be 


seen that 


the lift- 


drag 


function 


l remains 


relatively constant 


up to c/c 


= .58 

max 


Since 



the object is to make the span as small as possible this point 
is taken. The wing planform becomes, 
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c = .58 c = 3.62 ft = 43.44 in 

max 

c = 5.43 ft = 65.16 in 
r 

b = 2.02 ft = 24.24 in 
c = 3.90 ft = 46.80 in 
AR = 0.74 



5 . Design Example (For Canard Planform) 

From previous analysis and development, the following 
parameters have been defined: 

M = 3.0 

S = 1.09 ft 2 
c 

(t/c) c = 0.03 

If the leading edge sweep is set 5 degrees behind the sonic 



value and 


the taper 


ratio is 


set to zero. 


the following 


planform 


table results. 

a le 
X = 

AR b 


75.5° 

0 

c 

r 


c 


K W(b) 


B (W) 


2.5 


1.63 


1.34 


.89 


1.29 


.52 


2.0 


1.48 


1.47 


.98 


1.31 


.55 


1.5 


1.28 


1.70 


1.14 


1.35 


.59 


1.0 


1.04 


2.10 


1.40 


1.47 


.76 


.5 


.74 


2.94 


1.96 


1.51 


.90 


.25 


.52 


4.19 


2.79 


1.59 


1.02 
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Figure (7-8). Wing lift-drag function. 



AR 


% 


C °o 


AR tan 


S 


C Na (k > 


tan a le 


2.5 


.0069 


.0087 


9.67 


. 73 




2.0 


.0067 


.0085 


7.73 


.73 


2.7900 


1.5 


.0067 


.0085 


5.80 


.73 


3.6339 


1.0 


.0064 


.0082 


3.87 


.73 


2.5952 


.5 


.0061 


.0079 


1.93 


.73 


2.0880 


.25 


.0058 


.0076 


.97 


.73 


1.0124 



AR 


c/c 

max 


F 1 


F 

2 


F 


2.0 


.35 


1.0000 


.3628 


1.3628 


1.5 


.41 


1.000 


.3844 


1.3844 


1.0 


.50 


.9647 


.3901 


1.3557 


.5 


.70 


.9294 


.4849 


1.4143 


.25 


1.00 


.8941 


1.0000 


1.8941 



From Figure (7-9) it is seen that there is a minimum of 
the lift-drag function at c / c max = 0-5. The optimum canard is 
then , 

c = 0.5 c = 1.40 ft - 16.80 in 
max 

c = 2.09 ft = 25.08 in 
r 

b = 1.04 ft = 12.48 in 
AR = 0.99 

c = 2.01 ft = 24.12 in 



183 



1.9r 




-iT 



- 6 * 



- 8 * 



- r 



- 9* 



it — s* 



■ 



■ e* 



- 2 * 



- 1 



a 



to in 



m 



NOiiONDj ouaa-uii 



184 



Figure (7-9) . Canard lift-drag function. 



E. DESIGN EXAMPLE (REVISION OF ZERO LIFT-DRAG AND LIFT-CURVE 
SLOPE) 

Wing: c = 3.62 ft 

c f = .0014 

(C ) = .0750 

0 W 

Canard: c = 1.40 ft 

= .0016 

(C ) = .0166 

0 c 

Body: (C ) = .2292 

D 0 B 

The complete missile C including the interference factor 

0 

of 1.25, is now 0.4010. From the previous calculation of the 
baseline zero lift-drag coefficient, the drag has been reduced 
by 14.6 percent. 

From equation (6) in Chapter 6 the lift-curve slope is now 
C Na = 25.17/rad 

Figure (7-10) is the missile design to this point. 
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VIII. PROPULSION REQUIREMENTS 



The following discussion presents a method for preliminary 
sizing of a solid rocket motor for a boost-sustain trajectory 
of an air launched or a surface launched tactical missile. The 
analysis consists of sizing the booster from incremental velo- 
city considerations and sizing, the sustainer for the maximum 
range required at the operational altitude. The method assumes 
a constant acceleration boost and a constant altitude cruise. 



A. BOOSTER INITIAL SIZING 

Since the control system cannot respond properly while in 
the boost phase of flight, it becomes important to make the boost 
time as short as possible. The limiting factor for boost time 
is the maximum axial acceleration the airframe and components 
can withstand. This maximum acceleration is normally around 
30 g*s. The boost time then becomes. 




For constant acceleration, 





% 

a 



( 1 ) 



The incremental velocity during boost, AVg, can be derived 
from the equations of motion. 
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Figure (8-1) . Forces acting on the missile. 



If during the boost phase the thrust, drag and launch angle are 
considered constant, the axial acceleration of the missile is 
constant and may be written as: 

dV 

m -rr- = T - D - W sin y T 
dt ' L 

dV = g (— dt - g sin y dt 
W Li 



If the launch velocity is and the velocity at the end of 
boost is 




The vehicle weight in equation (2) is a function of time. If the 

propellant weight is given by W , and the propellant grain burn 

P 

is linear with time, the missile weight can be expressed as, 
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w = w. 



- W„(^) 

b 



L "p't 



where W = Launch weight. 

Li 

Equation (2) then becomes, 



av b = v 2 - Vj 



g(T - D) 



[ b t b dt 

J (W t -W t) sin y l b b 

n L JD p 



Integrating the first term yields. 



av b = 



g(T - D)t t 

w 



In 



% - 



w t 
Lb 



- 9 sin t b 



Since the empty weight is given by 



= W T - W 
E L p 



The incremental velocity can be written as 



g (T - D)t W 

av b = w £n wl - g sin y l b b 

P E 



(3) 



A first estimate of propellant weight can be made if the drag is 
assumed zero in equation (3) , and I = Tt^^ 



W L 

4V B = 9 1 Hi — - g sin y l 

£ 



(4) 



In the above expression the specific impulse, I , is charac- 

sp 

teristic of the propellant chosen and can be found from histori- 
cal data, or a specific propellant value may be used. 
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W L 

Equation (4) can then be solved for — and the first esti- 

W E 

mate of propellant weight, W , can be found as follows, 

P 



W E 

W p = W L (1 ' ST' 



With the propellant weight and specific impulse known, the 
total impulse and average thrust required for boost can be found. 



I 



T 

T 




W 

P 



With the thrust known and an average value of drag assumed, 
equation (3) can be iterated for an improved estimate of the pro- 
pellant weight needed. 

The final result is used to calculate the booster combustor 
volume and length. 



V B = W p/P p « p 



where n^ = volumetric packing factor 



L_ = 4V_ . ^ 2 
B B/ttD c 



The combustor diameter is limited by the missile maximum 

diameter and is a design choice. The propellant density, p , 

P 

and volumetric packing factor, n , are determined from histori- 
cal data or given for a specific propellant. 
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B. SUSTAINER INITIAL SIZING 



The sustainer thrust required to maintain cruise is the 

driving factor in the sustainer motor sizing. For initial sizing 

purposes, a level, constant velocity cruise segment is assumed. 

In this analysis the thrust required is equal to the aerodynamic 

drag developed by the missile. The performance requirements 

determine the maximum range, operating altitude, and velocity of 

the missile. From these requirements the sustainer burn time, 

t , can be determined, 
s 

, _ R max S B 

where S = distance traveled during boost. 

Since sustainer thrust is equal to drag, T = D, the total 
impulse required becomes. 







Dt 



s 



The sustainer nozzle is sized for the operational altitude 
of the missile. The thrust coefficient, C can be expressed as 



C 



Fd 



C d A 





Y+l 

Y-l 



1 




Y-l 

Y 



The constant, Cd, is the nozzle efficiency and has been 
determined from historical data to be 0.96. The nozzle half 
angle correction. A, is given by; 
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X = (1 + cos a)/2 



c*r 





where a is the nozzle half angle. A nominal value used in many 
designs is 15°. Larger values of a give significant non-axial 
flow components in the nozzle. Values less than 15° give large 
nozzle lengths and therefore contribute to excessive missile 
weights. The combustion chamber pressure is also a design choice. 
A lower usable value of P c is 200 psia for sustainers and 500 psi 
for boosters. The chamber pressure cannot exceed the maximum 
expected operating pressure of the missile (MEOP) . 

The thrust of the sustainer can be expressed as a function 
of the thrust coefficient, chamber pressure, P , and throat area, 

A . The throat area required to deliver the required sustainer 
thrust can be determined from 



T 

s 



'Fd 



The nozzle area ratio is a function of the pressure ratio 

P-/P , where P n is the ambient pressure if the nozzle expands 
u c u 

the flow fully. 





a±i 
v 2 



) M 




Y+l 

(Y-1) 



where 



P 



c 



1 + 






M 



Y 

Y-1 
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From the previous equation the exit area can be determined. 
From the throat and exit diameters, the nozzle length is found 
by 



L 

n 




2 tan a 



At this point in the design the nozzle exit area and length 
should be checked to see if they exceed any specifications on 
the missile; i.e., the nozzle exit area must be less than the 
base area of the missile. 

Finally, the delivered specific impulse is typical of the 
propellant chosen for the sustainer. The propellant weight is 
then, 



W = 
P 



(I T ) s 
I sp 



The sustainer combustor volume and length are then determined 
in the same manner as the booster. The equations concerning 
the sustainer nozzle also apply to the booster nozzle with the 
appropriate booster areas, pressures and thrust substituted* 

C. ROCKET MOTOR CHAMBER PRESSURE 

As indicated earlier the delivered thrust increases with 
increasing chamber pressure, P c ; however, the wall thickness 
in the rocket motor depends upon the expected operating pressure 
within the chamber. If the wall thickness is t inches, and the 
yield stress of the casing material is a , it can be shown that 
the thickness required of the motor casing is given by. 
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P r 

t = -SSS2 — (5) 

°y 

where r is the chamber radius, and P is the maximum chamber 

max 

pressure, which is taken to be 1.2 P c to allow for variations 
within the propellant. From equation (5) the weight of the motor 
casing can be determined if the casing is cylindrical and 
t < < r 



W = p V = 2nrl tp 
c c c c c 



Substituting for t 



W = 
c 



2.4irr SL p P 
c c c 



( 6 ) 



The specific impulse of the rocket motor is given by. 



= _1 = _i j 2g (— £y) R T 

sp g c g c ) c Y-l 0 



Y-l 
p 0 Y 
1 - 0 
c 



1/2 



(7) 



From equations (6) and (7) , it can be seen that a compromise 
must be made in selecting the chamber pressure; since increasing 
P c increases the specific impulse but also increases the casing 
weight. An attempt must be made to find an optimum operating 
pressure . 

If the rocket motor weight consists of the nozzle, propellant 
and casing, 



W„ = W„ + W_ + W_ 
M N P C 



( 8 ) 
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and the total impulse, which is constant is given by. 



I = Tt, = I W 
b sp p 



The optimum chamber pressure can be found by minimizing the 
weight and maximizing the impulse. 



- W M , W w 

d . _ d , M 

dp 1 I ; dp ( I W 
^c c sp p 



) = 0 



differentiating yields. 



w 

_i_ _a_ <:•«, 

r sp d Pc w p 



-L- _A_ . 

I dp V W 
sp C c p 



W M d ( 1/1 ) 

M sp 



w 

p 


dp 

c 


w„ 


d I 


M 


___ sp 


W I 


T d p c 



= 0 



P sp 



(9) 



From equation (8) 



W M W W 

= -£ + -N + 1 

WWW 
p p p 



(10) 



Substituting equation (10) into equation (9) , 

d I 

N £ _ M 

„ 2 dp W I 

W n sn 

P 



. dW 
1 c 

W dp 
p r c 



W dW 

c 



W 2 d ?c 



S£_ 



W I dp 
p sp ^c 



= 0 
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Substituting equations (12) and (13) into equation (11) 



2 . 4 p W 



W RT. 



£ EL E 2 (_£) 



Y-l 

Y 



p O g I ^ p p 

p y sp 



2 

Substituting for g c 1^ in the above equation yields , 



Y-l 


1 




2Y 


p n (1 -y)/Y 
(— ) -1 


P 




P c 


C 



2.4 



pc 

a 



0 



(14) 



If a propellant and casing material are chosen, equation (14) 

can be solved for the chamber pressures, p^, which minimizes 

W /I . One interesting feature of equation (14) is that there 

a y 

is a minimum yield strength to density ration — that will yield 

P c 

a usable chamber pressure, furthermore this optimization yields 
a very shallow minimum. Therefore, if the thrust or exit dia- 
meter requires an increase in p^ the penalty paid in additional 
casing weight will be small. 

D. TYPICAL PROPELLANTS 

From the preceding discussion it can be seen that some know- 
ledge of propellants to be used in the rocket motor is needed. 
The thrust developed by a rocket motor depends directly upon the 
pressure in the chamber. 



T = 




P 



C 
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Since 




I I 

sp sp 



dW 

c 

dp 



+ 



(W„ + W M ) 
c N 



I dl 
sp sp 



"sp 



d p. 



W M dI 
M sp 



sp 



dp. 



0 



From equation <8) 



W + W x = W„ - W 
c N M p 



Substituting for 



dW W dl 

_E R S P = o 

dp I dp 
c sp c 



From equation (6) 



dW 2. 4 nr l p 

c _ c c 

d P cr 



Since 



W = nr^i p 
P c K p 




2 . 4ttp W 
£ R 

P o 

p y 



from equation (3) 



dl 



££ = 



RT 



d P, 



g p I 
^c c c sp 



Y-l 
P o T 
(r 2 ) 



( 11 ) 



( 12 ) 



(13) 
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Because of this relationship, high pressures are needed for the 
boost phase, and relatively low pressure are needed for the sus- 
tain phase of flight. The pressure developed in a rocket motor 
chamber is a function of the burn rate and burn surface area of 
the propellant. In order to provide more burn area for the boost 
phase, the grain normally has an internal star perforation; while 
a sustain motor is normally an end burning grain of solid cross- 
section. The volumetric loading of a rocket motor is defined as 
the ratio of the propellant volume to the rocket motor chamber 
volume . 



Due to erosive burning effects the volumetric loading of a booster 
is normally limited to less than 0.9. For efficient packing 
the ratio varies from 0.7 to 0.9. The volumetric loading of an 
end burning sustainer engine is normally 1.0. The range of 
propellant characteristics is shown in Table (8-1). 

TABLE (8-1). Typical Propellant Properties [14,15]. 



_ Grain volume 
p Chamber volume 



Sustainer 



Booster 



I „ (sec) 
sp 



180-210 



210-260 



p (lbm/in^) 



P 



059- . 062 



062-.065 



Y 



1.24 - 1.27 



1.22 - 1.26 



n 



1.00 



0.85 



P 
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E. DESIGN EXAMPLE (BOOSTER) 

From the previous analysis the launch speed is M = 1.5, 
the cruise speed is M = 3.0 at an altitude of 10,000 feet, 
boost incremental velocity is then. 



AVg = 1616.10 ft/sec 

If the maximum acceleration during boost is a constant 30 gs 

(assume sea level g ), the time of burn is, 

c 



t 



b 



iV B /a 



max 



1616.10 ft/sec 
30(32.2 ft/sec 2 ) 



t^ = 1.67 sec 



From equation (4) with y = 0 

1j 

. «L iV B 
n W E " 3 



Assume I = 250 sec from Table (I) 
sp 



W L 

l rF- = 0.2008 

n W E 

W F 

W = W_ (1 - A 
p L Wg 



For 



W T = 794.83 
Li 



W = 144.57 abm 
P 



and 

The 
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The thrust provided is 



I W 

(T). = — 2- = 21,642 lb, 



Since the missile flies essentially at zero angle of 
attack during boost, the drag at the end of boost is. 



D = (C ) q S 
u 0 



ref 



D = 2006 lb J 



Therefore if an average drag of 1000 lbs is assumed, equation (3) 

can be used to find a new W./W^, 

L t, 



w. 



AV_ W 
B 



'n W E g(T - D)^ 




.2270 



W = 161.40 lb 



p m 

This gives a thrust of 24,162 lbs. 
equation (4) gives, 



One more iteration of 



W = 159.46 l & 
p m 

(T) = 23,872 lb £ 

s f 
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Continued iteration of equation (4) does not change the pro- 
pellant weight. 



From Table (I); p = .062 Jtb /in 3 

p m 

n = .85 
P 

The booster volume then becomes 



II 

> 


3025.81 in 3 


Allowing one-half 


inch for structure, the booster length becomes. 


l b = 


4 V B 

= 42.69 inches 

ttD 

c 



F. DESIGN EXAMPLE (SUSTAINER) 

The maximum- range of the missile was determined to be 35 nmiles 

R = 212,800 ft 
max 

The distance traveled during boost is 



S B = 


V M S 


V = 


V 1 + V 2 

— - = 2424.15 ft/sec 


S B = 


4048.33 ft 



Therefore, at a cruise Mach number of 3.0, the sustainer time 
of burn is 

t = 64.59 sec. 
s 
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The cruise drag consists of zero-lift drag and induced drag. 



The lift coefficient of the missile minus the booster is; 



= 



W 



L q S 



ref 



C T = 0.1270 
Li 



From Chapter 6 



C I> C D Q + C L /C Nct 



C D = .4016 



The cruise drag then becomes 



D = 2009 £b f 



The total impulse required is then; 

(I_) = Dt = 129,761 Zb £ - sec 

T s s r 



From Table (I) the specific impulse is 210 sec; therefore, 
the weight of propellant needed for cruise is 



(I T } s 
(W ) = - T ^ 

p 3 I sp 



= 618 Zb 



m 



The volume required is 

V = 9967.7 in 3 
s 
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The sustainer length becomes, 



4V 



L 



s 



^ = 140.62 inches 



c 



The total rocket motor length is, 




L 



RM 



183.31 in 



The rocket motor alone exeeds the total length of the base- 
line missile. This will lead to large length to diameter ratios 
for the entire missile; therefore, the rocket motor length must 
be decreased. That is the subject of the next section. 

G. REDUCING ROCKET MOTOR LENGTH 

If the rocket motor length of the last section leads to 
excessive length to diameter ratios for the missile, the motor 
length must be decreased. This may be accomplished in one of 
two ways: (1) If the mission profile has a long cruise segment, 

the cruise altitude may be increased thereby decreasing the drag 
and the total impulse required. (2) If the cruise altitude 
cannot be varied, the missile diameter must be resized for pro- 
pulsion considerations. In this design example option (2) will 
be selected. 

The missile length can be expressed as 




( 15 ) 
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L n = Nose length 
L_ = Guidance section length 
L c = Control section length 
L WH = War ^ ea< ^ section length 
L RM = Roc ^ et motor length 



The nose and warhead sections have been previously defined. 
Since the propellant weight for the range requirement is known, 
the diameter necessary can be determined if a maximum missile 
length is specified. The drag will increase slightly due to the 
increased diameter but the total impulse, and therefore the pro- 
pellant weight, will change only slightly. 

The rocket motor length can be expressed in terms of the 
diameter by 



4W 



RM 



V° 



Equation (15) then becomes 



(L) = L n + L g + L o + L wh + 



4Wp 



PpTTD 



The guidance and the control sections are each normally about 

10 percent of the missile length. 

L L 4W 

(L) = (— )D + . 1L + . 1L + (-—) D + — E-y 

D D p ttD 

P 



n 



WH 



- . 8LD + 



4W 

_ P., = 

P 7T 

P 



(16) 
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As mentioned in Chapter 2 the missile is normally designed for 

an existing launcher. This launcher will have a maximum length 

that it can accomodate, L . If this value is substituted into 

rna.x 

equation (16) , the required diameter can be found. 

1 . Design Example (Resizing for Propulsion) 

From previous analysis; L^/D = 2.23 

L WH //D 2,50 

The propellant weight from before was, W = 777.46 lbs. If 

P 

the maximum length of the missile is taken as 210 inches, equa- 
tion (16) becomes, 

4.73D 3 - 168D 2 + 15,467.08 = 0.0 

The required diameter is then, 

(D) = 11.75 inches 

req 

Allowing 0.25 in for structure, the missile diameter becomes 
12 inches. From the equations for the rocket motor lengths. 



= 26.62 in 
L = 90.47 in 

b 

L = 117.09 in 
RM 



H. DESIGN EXAMPLE (CHAMBER PRESSURE) 

From equation (14) , the strength to density ratio can be 

solved for in terms of p . 

^c 
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a ' 

_X 



. 2y 



(Y-l) 



min 



1-Y 

P o Y 
( P £) 



-1 



P <^±) 
*c p ' 



If the minimum usable chamber pressure is = 500 psi, and 
using the propellant previously selected. 



Y = 1.24 
p = .062 

P 

P 0 = 10.11 



the minimum strength to density required can be found. 

a 

(-£) = 225,542.87 in 

p c 

min 

Inconnel 718 is a common metal used in combustion chambers 
and has the following properties [4]: 



a = 180,000 psi 

y * 

p = .2662 b /in 3 
c m 



-2- = 676,183 > (-£) 

) p 

c c 



mm 



The qtimum chamber pressure is then from equation (14) 



P c = 1132.55 psi 



The rocket motor chamber wall thickness is given by equation (5) 
as 



206 



t 



r 



P. 



max 



t = 0.05 inches 



The wall weight is given by equation (2) 




W = 59.5 Jib 
c 



m 



I. DESIGN EXAMPLE (SIZING THE NOZZLE) 

For initial analysis it is assumed that the flow is fully 
expanded, in a nozzle with a half angle of 15 degrees. With 
the chamber pressure of 1132.55 psi the thrust coefficient at 
altitude is 




C 



F 



2.3023 



d 



The throat area and diameter then become 



A 



T 



9.1552 in 



2 



t C 



P 



c 



F 



d 



d 



t 



3.4142 in 
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From isentropic tables for Pg/P c = .0089 



d 




1/2 



d 



. 3342 



The exit diameter becomes 



d = 10.22 in 



e 



The nozzle length for a 15 degree half angle is. 



l = 12.70 in 
n 



The nozzle is shown in Figure (8-2) . Figure (8-3) is the design 
with resized diameter and length for propulsion and the optimum 
wings determined in Chapter 7. 

J. TAIL SIZING 

As mentioned in Chapter 7 the tail (canard) sizing is accom- 
plished for zero static stability at launch. Now that the missile 
length has been fixed, the canard can be sized. Figure (8-3) 
is a drawing of the design to this point. The total lifting sur- 
face required is 



S TT + S = 6.57 ft 
W c 



2 



if 




( 16 ) 



From equation (5) in Chapter 7, for C = 0 
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c 

•H 

<N 

CM 

O 



<D 

"3 
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Figure (8-2) . Nozzle geometry. 



40 in 



•f* 







210 



0 



(C Na> + (C Na> x w s" 2 , (K W, k . + K b, w > 
N W ref (B) (W) 



(17) 



+ (C.. ) x — 
Net ' C S - 
C ref 



(K + K ) 
U (B) B (C) 



Equations (16) and (17) above can be solved for S TT and k, 

w 



S w = 4.70 ft' 



k = 0.40 



Therefore, S =1.87 
c 



If the same aspect ratio as found in Chapter 7 is used, the wing 
and canard planforms become. 



: b = 1.86 

w 


ft 


Canard: b = 

c 


1.36 


ft 


(c ) = 5.05 

r w 


ft 


( Vc ■ 


2.75 


ft 


(c) = 3.60 

w 


ft 


(c) c = 


2.63 


ft 
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IX. NONLINEAR AERODYNAMICS AND AEROl 



As mentioned in Chpater 5 the analysis thus far has been for 
one speed and small angles of attack. The aerodynamics change 
dramatically from subsonic to supersonic flow and with increasing 
angle of attack. This chapter presents the methods used primarily 
by USAF Stability and Control DATCOM and covers all configura- 
tions and flight speeds. As will be shown, the method becomes 
very involved and therefore lends itself to programming on a 
digital computer. 

The lifting characteristics of both wings and bodies become 
nonlinear as the missile angle of attack increases above 10°. 

Dp to this point in the thesis development only the linear 
contribution has been considered. For large angles of attack, 
the nonlinear term, which is due to flow separation, must be 
considered. The relative effect of these terms on the normal 
force and pitching moment coefficients is shown in Figure (9-1) . 



A. VISCOUS CROSS-FLOW 



As can be seen in Figure (9-1) , at large angles of attack 

the forces acting are primarily nonlinear. The nonlinear term 

is normally described through the use of the viscous cross- flow 

coefficient, C. 

d 




v 



q S 
M n p 



Sp = planform area. 
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j 
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Figure (9-1) . Relative contributions of linear and nonlinear 

terms [ 19 ] . 



For an infinitely long circular cylinder, the viscous cross- 
flow force per unit length is 

f = 2r C , q n (1) 

v d c n 

The viscous cross-flow coefficient is determined experimentally 
and is a function of normal Mach number and Reynolds number. 
Figure (9-2) gives the cross-flow drag coefficient as a function 
of cross-flow Mach number. 
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Figure (9-2). Cross-flow drag coefficient [18]. 



M = M,. sin a 
n M 



R = R sm a 
e e 

n 

V = V M sin a 
n M 



Since in real flow there is spillage around the ends of a finite 
length cylinder, the drag is less than that for an infinite, 2-D 
cylinder. The cross-flow drag proportionality constant, n, is 
the ratio of the drag coefficient of a finite cylinder to the 
drag coefficient of an infinite cylinder. 



n = 



(C. ) 
d c 2-D 



The proportionality constant is also determined experimentally 
and is given in Figure (9-3) . Notice that the cross-flow drag 
proportionality constant approaches one as the length to dia- 
meter ratio becomes large, or the 2-D situation is approached. 

The viscous cross-flow force per unit length for a finite 
cylinder then becomes 

pv 2 

f v - 2rnc d — 2 s - t2) 

c 



Since 





. 2 
sin a 



. 2 

q sin a 



equation (2) becomes. 



f 



v 



2n V 



. 2 

r sm a 
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The viscous contribution to the normal force acting on the cylinder 
then becomes 



(C N> 



V 



q s 



ref 



2nC . sin a /.£ 
d c 

S _ 
ref 



rdx 



For a uniform cylinder, 



n C d S f 
V c ref 



P • 2 
£ — sm a 



If this force is added to the term predicted by slender body 

potential theory, the total normal force acting on the body can 

be found. Reference [19] gives this total force as; 

S, -> S 

,_ N b . 2 a . p . 2 

(C N ) = sin a cos ^ + n C . — ^ ■ sin a 

B S ref c b ref 



Similar terms can be developed for the wing and tail since 
they also exhibit nonlinear behavior at high angles of attack. 
Care must be taken to separate the lift of Chapter 5 into the 
lift acting on the wing panel and the additional lift on the body 
before the cross-flow term is applied. The nonlinear cross-flow 
component of wing lift is caused by flow separation and the for- 
mation of vortices on the upper surface of the wing. This 
viscous term is given by, C 



'N. 



W(V) 



. 2 w 

: = C, sin a ^ 

N T7/TT . d S , 

W(V) n ref 



Where C^. is the cross-flow drag coefficient for a flat plate and 
n 

is given in Figure (9-4) . 
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B. TOTAL MISSILE LIFT 

The remainder of this chapter is a summary of references [17] 
and [18]. With the addition of the force due to viscous cross- 
flow, the total lift on the body becomes 




(( V 



+ C 



B 



N 



+ C 



B (W) 



N 



+ C 



B (T) -a 



N 



B (T) -6 



where , 

(C M ) = linear and nonlinear lift force on the body 
N B 

C . = additional lift on the body due to the wing 

N B(W) 

C. T = additional lift on the body due to the tail 

N B(T)-a 

(due to angle of attack) 

C„ = additional lift on the body due to the tail 

N B(T)-<5 

(due to control deflection) 



The additional lift on the body in the presence of the wing and 
tail can be found as described in Chapter 6 . 

The wing lift is now composed of a linear and a nonlinear 
component and becomes 



C„ = C M + C XT 

N W N W(B) N W(V) 

S 2 

C N T7 = (B) (C N**\ 7 S 7 sina + C d Sin a S 

W W ref n 



The tail lift is found in a similar manner only with the 
additional term for any control surface deflection, 6. 

N T N T(B)-a N T(B)-S N T(V) 
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Figure (9-4). Drag coefficient for a flat plate normal to flow [18]. 



Notice that the linear terms in the above equations are multiplied 

by sina instead of a as in Chapter 6. The linear theory from 

which the Chapter 6 equations were derived assumes small angles 

of attack and therefore sina^a . This is not true at higher 

angles of attack and is therefore included here. 

The total normal force then becomes , C„ 

N 



'N 



= C 



N 



+ C 



B 



N, 



+ C. 



W 



N„ 



+ C 



L i 



The last term, C T , is the lift-loss due to downwash and is given 

Li . 

1 

by 



' La w La T 



W 



(B) 



W 



(B) 



(C T „,) (C T ^) |~K rT sina + k rT sinSj i (b - r) m S 



T w 



'L. 

l 



2n AR m (f - r ) S - 
T w w ref 



' i = interference factor 

f = votex location 
w 

the above terms are found by the methods of reference [13]. 

C. DRAG CHARACTERISTICS OF A MISSILE 

The total drag acting on a missile is the sum of the zero- 
lift drag and the induced drag (due to angle of attack and/or 
control surface deflection) . The zero lift drag of bodies and 
wings is highly dependent on the missile speed. 

Three speed regimes are normally considered. 
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Subsonic 0 < M„, < .8 

- M 

Transonic .8 < M w < 1.2 

- M - 

Supersonic , • M > 1.2 

Because of the empirical nature of the formulas for zero-lift 
drag, the procedures followed are essentially those of DATCOM [16] 
1 . Zero Lift Drag 
a. Subsonic 

At subsonic speeds the zero-lift drag consists of 
skin friction (incompressible) and pressure drag. The pressure 
drag at subsonic speeds is usually small compared to the drag 
due to skin friction. The entire zero-lift drag of a missile 
at subsonic speeds is given by; 



= c: + 

D 0 °0B D 0W d ot 



Where the components are the body, wing and tail contributions 



(1) Body drag, C 



D 



OB 



The body zero-lift drag is 



given by; 



'D 



= 1.02 r 



OB 



-B 



1 + - 1 ^ + — 
f 3/2 + f 3 



(S Wet ) 



B 



+ C, 



ref 



B 



where C is found by the methods of Chapter 5. 

A B 

The first term is the skin friction contribution and the 
next two terms are the pressure contributions. f is the body 
fineness ratio and is given by; 
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(2) Wing drag, C . The wing zero-lift drag is 



OW 



given by. 



C = C, [l + 2(|) + 100 (|) 
D 0W f w L c c J 



4 1 (S Wet ) 



W 



ref 



^ S Wet^ ; '' s ent ^- re w ®tted area of the wings. The thickness 



W 



to chord ratio is given by (t/c) 



(3) Tail drag , C . The tail zero-lift drag is 

°0T 



given by. 



0T 



a C- fl + 2 (-) + 100 (-) 

f T L C c J 



41 (S Wet\ 



ref 



The above analysis assumes fully turbulent flow along all sur- 
faces. , C and C_ are the local average skin friction 

£ B r W r T 

coefficients based on the local Reynolds number, . 

The reference lengths are the body length and 
wing/tail mean aerodynamic chords. The skin friction coefficient 
is then given by 



C 



f 



4.55 

UogiQ R e > 



2.58 



b. Transonic 

The total transonic zero-lift drag is composed of 
skin-friction drag, transonic wave drag, pressure drag and 
base pressure drag. The compressible skin friction drag is 
found by using a correction factor on the incompressible co- 
efficient, , found for subsonic flow. The compressibility 
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correction is found by the methods of Chapter 6. The skin 
friction drag then becomes 



(S Wet ) 



C D 1,02 C f s ~ 

fB c ref 



B 



The subsonic pressure drag is the same as before 



C = 1.02 C f 
pB r B 






5 + 7 

3/2 + f . 



(S Wet } 



B 



ref 



The above equation applies for Mach numbers in the range of 

.8 to 1.0. The pressure drag then decreases linearly to zero 

at a Mach number of 1.2. The transonic wave drag is deter- 

rB 

mined using Figure (9-5) , which is a function of nose fineness 
ratio L N/D * 

The total transonic zero-lift drag of the body then becomes 



OB 



= C + C + C 
D. DpB D 



N 

VB S ref 



+ C, 



B 



where S N is the cross sectional area at the nose-body junction. 

The total transonic zero-lift drag of the aerodynamic surfaces 
is composed of the skin friction drag and a drag increment, 

AC d , which is the transonic wave drag. 

Experimental results show little increase in the viscous 
drag of aerodynamic surfaces from the subsonic to the transonic 
regimes; therefore the skin friction and pressure drag is the 
same as for subsonic flow and is given by 
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c 



D 



OW 




1 + 2<i> 



+ 100 <-) 
C 



S 

S 



Wet 

ref 



The wave drag of the aerodynamic surfaces if found from 
Figure (9-6) and is a function of AR and M. For swept wings 
the Mach number used in Figure (9-6) is given by. 



M ' = M [cos A c/4 ] 



1/2 



The transonic wave drag increment is then given by 



2.5 



AC. 



= AC. 



[cos A c/4] 



OW 



OW 



W 

’ref 



where AC is obtained from Figure (9-6) . 

D OW 

drag is found in the same manner as the wing, 
drag (transonic) is then given by; 



The tail zero-lift 
The total zero-lift 



= C. 



+ C. 



+ AC. 



+ C. 



+ AC. 



OB 



OW 



OW 



OT 



0T 



c. Supersonic 

The supersonic zero-lift drag of a missile is deter- 
mined empirically by assuming a parabolic variation of C with 

D 0 

Mach number between 1.2 and 3.0. The resulting equation is for 
the entire missile and is given by 

C = C ' - 1.7209 (C " - C ') + 1.5708 (C ” - C ') /m 
0 0 D 0 D 0 U 0 U 0 
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l n/d 



Figure (9-5) . Transonic wave drag for ogival and blunted 

conical forebodies [18]. 
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Figure (9-6). Transonic C for unswept wings [18]. 
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Figure (9-7). Ratio of C D for various L/D to C D hemispherical [18]. 



where 



C 



I 




(C^ ) 



M 



1.2 



C 



D 



If 

0 



(c 




M 



3.0 



C ' is determined by using the methods for transonic flow. 

D 0 

The magnitude of the supersonic wave drag is highly 

dependent upon the nose shape. For this reason two methods 

are used to determine C " 

0 

Method I : This method is for blunted ogives, pointed ogives and 
blunted cones. In this method 



For Vd * * 5 

For L N/D i 8 • ' 0 





(C dJ 



M 



+ AC r 



.8 



0W 



+ AC t 



+ C, 



0T 



B 



The above values are determined from the methods of transonic 
flow. 



For L 
C 



N/D 



It 




.5 to 8.0 



K 



1 



C 



f 




in the above equation is derived empirically and is given in 
Figure (9-7) . 

Method II : This method is for pointed conical noses. In this 

method C D " is determined by 
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'D, 



= (c dJ 



+ AC 



D 



+ AC. 



+ (c r 



+ c. 



M = .8 



OW 



OT 



VN 



M = 3 



B 



The first three terms in the above equation are found by transonic 

flow methods. The forebody wave drag, C , for M = 3.0 is found 

D VN 

using Figure (9-8). The nose semi-vertex angle, 8 , is the same 
as a in Chapter 6. In all flow regimes the base pressure drag 
is found as in Chapter 6. 

2 . Induced Drag 

The induced drag due to angle of attack depends upon the 
flight regime the missile is in. For M < .85 and AR > 3.0 the 
induced drag is, 

(C L )2 

C = — 

D . TTARe 

l 



where e is the Oswald efficiency factor » 0.7. 

The tactical missile normally has an aspect ratio of less than 
3.0. For all flight speeds with AR < 3.0 the induced drag is 



C_ = C T tan a 
D • Li 

1 

The component induced drag coefficients are found in the same 
manner as above using component lift coefficients. 

3 . Total Drag 

From the zero-lift drag and induced drag the missile 
total drag becomes , 



C 



D 
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Figure (9-8). Wave drag of pointed conical nose [18], 



D. PITCHING MOMENT CHARACTERISTICS 



The pitching moment of a missile is equal to the sum of 
the pitching moments due to lift and drag forces acting on the 
body, wings and tails. If only small angles of attack are of 
interest, the pitching moment is due primarily to the lift 
forces. The methods presented here are for all angles of attack. 

1 . Body 

The body alone pitching moment about its- center of gravity 
is given from slender body theory and viscous cross-flow theory 
as 





- S B (L - X cG ) 
S ref L ref 



. 2 

sm a cos 





S x _ - x ~ 

( _c| — ) sin^a 

ref ref 



S = Planform area 
P 

Xp = Centroid of planform. area 

V = Body volume 
B 

S = Base area 
B 



With C^ given the center of pressure of the body becomes 



B 



x 



(x = > ■ 
P B 



cG 



'M. 



B 



J ref 



'N 



J ref 



B. 



2 . Wing (Fixed Surface) 

The center of pressure location for the wing must be 
found before the pitching moment can be specified. The 
center of pressure of the various lift components are found 
by the method of Pitts, Nielsen and Fhatiari [13]. The lift 
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components are the same as in the lift section. The center 

of pressure of the lift on the wing in the presence of the 

body, (x) , is found from Figure (9-10) if the flow is 

W (B) 

subsonic and Figure (9-11) if the flow is supersonic. 

The center of pressure of the additional lift on the body 

in the presence of the wing, (x) , is found using Figure 

B (W) 

(9-12) if the flow is subsonic. If the flow is supersonic 
the center of pressure is found from Figure (9-13) or Figure 
(9-14) depending on the parameters; 



BAR ( 1 + A) (1 + 






A = taper ratio 
m = c° + 4 le 



If the center of pressure reference is moved to the nose, the 
following expressions result. 



(x 

(x 



cp J W 
cp' 1 B 
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Figure (9-10) . Wing center of pressure 

(subsonic) [ 18 ] . 
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Figure (9-11) . 
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Figure (9-12). Body center of pressure [18]. 
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Figure (9-13) . Body center of pressure [18] . 
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Figure (9-14) . Body center of pressure [18] . 
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From the above figure the entire pitching moment about 
the nose of the missile due to the forces acting on the wing is. 






C = [ (C + C + C )cosa + AC sina] ( _ 

”w L W(B) L W(V) L i D w-a L 



(B) 



ref 



(x ) 
cp 



+ [C cosa] ( 

L B(W) L ref 



B (W) 



The viscous lift, C and lift-loss due to downwash, 

L W(V) 

C , are not shown in the figure but are the same forces as 
L . 

1 

found in the lift section. These forces are assumed to be 

acting through the center of pressure, (x ) 

Cp W(B) 

An average center of pressure of the wing due to all forces 
acting on it can now be found as 

C M L « 

. . _ “w ref 

X cp T , Tel + C T ] cosa + AC sina 

W L t , L . DrT 

W i u W-a 



where 

L W L W(B) L B(W) L W(V) 



3 . Tail (Control Surface) 

The tail pitching moment about the nose of the missile 
is found in the same way except now a control surface deflec- 
tion must be included. The equations now become 
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(V) 




[ (C T + Cl 

T(B) T 



+ C ) cosa + C + 

L i L T(B)-S 



(x ) 

(AC + AC ) sina] ( — ) 

°T-a D T-6 ref 

(x ) 

+ [C cosa + C ] ( — 2JZL) 

^B (T) B (T) -<5 ref 

C L _ 

M ref 

( x ) = ± 

cp (C T + C )cosa + C + (AC + AC ) sina 

T ^T-a L i L T-5 D T-a D T-<5 



where , 



C = C + C + C 

L T-a L T(b) L b(T) L T(v) 



C = C -f C 

L T-S L T (B) -5 L B(T)-5 



The wing and tail pitching moments above were found about 
the nose of the missile. Transferring the pitching moments 
about the center of gravity the complete missile pitching 
moment becomes 



C = C 
M M 



+ C 



B 



M, 



W 



'x - (x ) 

C G cp 


-u r 


x - (x ) 

cG cp T 


(X ) 

P W 




(X ) 

rp 



E. AER01 DESCRIPTION 

It can be seen from the preceding description of the aero- 
dynamics of a missile, that the process of obtaining a complete 
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description of the aerodynamic coefficients becomes an involved 
task. This was the justification for initially using linear 
aerodynamics. To complete an analysis a fast method of pre- 
dicting the static aerodynamic characteristics is needed. 

The method needs to include both linear and nonlinear contri- 
butions as well as interference factors and must be applicable 
to all speed regimes. This analysis lends itself to program- 
ming on a digital computer. AEROl is a modification of the 
program AEROCF which was developed at the Naval Air Development 
Center by Mr. F . A. Kuster, Jr. [17] . The program is essentially 
as written except for calculation of planform areas and centroids. 
The program was also modified for use on the Naval Postgraduate 
School CP/CMS system. 

AEROl consists of a main program and three subroutines. The 
inputs to the program are the geometric characteristics of the 
missile, the flight conditions, engine and inlet type and the 
protuberance drag. The output consists of the aerodynamic co- 
efficients for lift, drag and moment. The component contribution 
to these coefficients are also given. The component and overall 
center of pressure are also determined. 

Subroutine GEOSUB ; This subroutine calculates the missile 
wetted area and the Reynolds number per foot based on the flight 
altitude . 

Subroutine CLASUB ; This subroutine calculates the aerodynamic 
surface lift-curve slopes. 
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Subroutine CATSUB; This subroutine calculates center of 



pressure locations, cross-flow drag coefficients, and interference 
factors . 

Main Program ; The main program assumes the control surface 
is the tail. This is regardless of the method of control (Wing, 
Tail, Canard) . Because of this, care must be taken to input the 
right data for the tail. For instance, if the missile is a wing 
control missile, the wing data is input as the tail and the 
tail data as the wing. Figure (9-15) and Table (9-1) give a 
complete listing of the input data. Table (9-II) is a list of 
the output data. Appendix A is a listing of the program as modi- 
fied for use on the Naval Postgraduate School IBM 360 computer. 

1 . Verification of AER01 

Before using the program an attempt was made to verify 
its prediction techniques and find any pitfalls in its use. To 
accomplish this the program was run for various configurations 
for which experimental data were available, and the results were 
compared. The comparisons are shown in Figures (9-16) to (9-25) 
from references [19] - [22] which are NASA technical notes and 
memoranda which report results of wind tunnel tests for various 
body-wing- tail combinations. In Figures (9-16) to (9-25) the 
solid lines are AEROl predictions, 
a. NASA TN D-6996 

This technical note presents a method of predicting 
aerodynamic characteristics for bodies alone at angles of attack 
from 0 to 180 degrees. Several nose-body combinations are given. 
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TABLE (9-1) 



Input Data 



Name 


Format 


Meaning 


ICSC 


12 


Type of control (wing, 
tail, canard) 


INOSE 


12 


Nose shape (ellipsoidal, 
cone, ogive) 


IDT 


12 


Number of control surface 
deflections 


IM 


12 


Numer of Mach numbers 


IAL 


12 


Number of angles of attack 


NBODY 


12 


Number of configurations 


ISWPW 


12 


Wing shape (delta , nondelta) 


IAFBW 


12 


Missile body after wing 


IWEPW 


12 


Leading edge sweep indicator 


NWING 


12 


Number of wings 


ISWPT 


12 


Tail shape (delta , nondelta) 


IAFBT 


12 


Missile body after tail 


ISWEPT 


12 


Leading edge sweep indicator 


NTAIL 


12 


Number of tails 


XLAMW 


F10 . 5 


Wing taper ratio 


CLAMW 


F10. 5 


Wing leading edge sweep 


BW 


F10 .5 


Wing span 


CROOTW 


F10. 5 


Wing root chard 


SW 


F10. 5 


Wing exposed area 


XMACW 


F10.5 


Wing mean geometric chord 


XWING 


F10.5 


Distance to wing leading 
edge 
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Variable 


Name 


Format 


Meaning 


22 


TOVCW 


F10 .5 


Wing thickness to chord ratio 


23 


XLAMT 


F10 .5 


Tail taper ratio 


24 


CLAMT 


F10 .5 


Tail leading edge sweep 


25 


BT 


F10 .5 


Tail span 


26 


CROOTT 


F10 .5 


Tail root chord 


27 


ST 


F10 .5 


Tail exposed area 


28 


XMACT 


F10. 5 


Tail mean geometric chord 


29 


XTAIL 


F10 .5 


Distance to tail leading 
edge 


30 


TOVCT 


F10 .5 


Tail thickness to chord ratio 


31 


HT 


F10 .5 


Altitude 


32 


D 


F10 .5 


body diameter 


33 


XL 


F10 .5 


Body length 


34 


XLNOSE 


F10. 5 


Nose length 


35 


XCG 


F10.5 


Center of gravity location 


36 


AREA 


F10. 5 


Reference area 


37 


XREF 


F10. 5 


Reference length 


38 


ENGINE 


F10.5 


Engine code 


39 


INLET 


F10.5 


Inlet code 


40 


BETA 


FlO .5 


Boattail angle 


41 


DBASE 


F10.5 


Base diameter 


42 


DJET 


FlO. 5 


Nozzle exit diameter 


43 


XLABOD 


FlO. 5 


Boattail length 


44 


CDPROT 


FlO . 5 


Proturberance drag 
coefficient 
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TABLE (9— II) 



Output 


Data 


Variable Name 


Meaning 


AL 


Angle of attack 


CLTOT 


Total coefficient of lift 


CDTOT 


Total coefficient of drag 


CLWP 


Wing panel coefficient of lift 


CLBW 


Additional lift on body due to 
wing coefficient 


CLTP 


Tail panel coefficient of lift 


CLBT 


Additional lift on body due to 
tail coefficient 


CLB 


Body alone lift coefficient 


CDI 


Induced drag coefficient 


CNWP 


Wing panel normal force coefficient 


CNTP 


Tail panel normal force coefficient 


CLTD 


Coefficient of lift due to tail 
deflection 


CDTD 


Coefficient of drag due to tail 
deflection 


CN 


Total normal force coefficient 


CA 


Total axial force coefficient 


XCPW 


Wing center of pressure 


XCPT 


Tail center of pressure 


XCP 


Total missile center of pressure 


CM 


Total pitching moment about center 
of gravity 
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Figures (9-16) to (9-18) compare the coefficients predicted by 
AER01 with those obtained for body number 9 in the NASA report. 

The normal force coefficient is predicted well throughout the 
range. The moment and axial force trends are predicted by the 
program but large errors exist throughout the range of angles 
of attack. 

b. NASA TM X-2367 

This technical memorandum investigates the aero- 
dynamic characteristics of various cruciform body-wing combina- 
tions. The coefficients for this configuration agree very well 
with experimental values up to 10 degrees angle of attack. Beyond 
this value the lift and moment coefficients predicted by AEROl 
exceed the experimental values by as much as 10 percent. Al- 
though the exact cause of this error was not investigated, it may 
be partially explained by the nose shape of the body. The nose 
is a combination ogive and cone. For purposes of AEROl it was 
assumed an ogive. The greater presented area of the ogive would 
contribute to a higher and C m through the cross-flow terms in 
these coefficients. Figures (9-19) through (9-21) compare the 
predicted with the experimentally determined coefficients for 
the wing-tail configuration of this reference at M = .9. 

c. NASA TM X-2780 and NASA TM X-2289 

These technical memoranda investigate the aerodynamics 
of a delta wing missile using tail control and a tail-less cruci- 
form missile. As shown in Figures (9-22) to (9-25) there is 
excellent agreement between the experimental values of the aero- 
dynamic coefficients and those predicted by AEROl. 
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Figure (9-16). Body alone normal force. 
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Figure (9-17) . Body alone moment. 
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Figure (9-18). Body alone axial force. 
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Figure ( 9 - 19 ). Tail-body lift force. 
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Figure (9-20). Tail-body drag force. 
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Figure (9-21) . Tail-body moment. 
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Figure (9-22) . Wing-body normal force. 
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Figure (9-23) . Wing-body moment. 
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Figure (9-24). Complete missile lift force. 
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Figure (9-25). Complete missile drag. 



F. COMPONENT WEIGHTS 

As with the gross weight, the component weights can be 
determined through the use of parametric regression equations. 
If the dimensions of the components are known, the following 
formulas can be used to determine the component weights [7]. 
Aero Surface Weight 

The weight of one wing panel is given by, 

W AS = 6.77483 (E^) 1 * 02 (AR^)* 56 (1) 

2 

E Exposed are of one panel, ft 
AR as Aspect ratio of one panel 

Body Structure Weight 

W BS * ' 0604 'W 64 (D BS )1 ' 77 t2 > 

L Length of body to be covered by 
BS 

structure. This does not include 

the rocket motor unless a separate 

structure surrounds the motor 

casing. (inches) 

Diameter of body structure 

°BS 

(inches) 

Internal Systems Weight 

»IS - - 00485 <V 74 H-is ’ 1 ' 00 (D IS>' 42 (3) 

W_ Gross weight of the missile 
G 
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L IS Length of subsystem (inches) 
D^g Diameter of missile (inches) 



In the above equation, the internal system includes guidance, 
avionics and control. 

1. Design Example (Component Sizing) 

If the guidance and the control sections are kept at 
10 percent of the total length, and the remaining components 
are as sized previously, the design is as shown in Figure (9- 
Since these components are considered internal systems, their 
weights can be determined from equation (3). 



W = .00485 ( 1809 ) " 79 (.1L) (12.5) * 42 

cont 



W ^ = W . , = 83.64 lb 
cont guid 



The wing and canard weights are found using equation (1) and 
Figure (9-26) . 




Figure (9-26) . Aero surface weight. 
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The wing weight is four times the weight of one panel. 

W T7 = 4 3.42 lb 
W 

Similarly the canard weight is 

W - 9.88 lb 
c 

If the entire body is covered by structure, the body structure 

weight becomes 

W = .0604 (232) * 64 (12. 5) 1 * 77 
st 

W = 172.38 lb 
st 

The following component table can now be constructed. 



Component 


Length (in) 


Weight (lbs) 


Center of Gra 


Control 


23.2 


83.64 


39.48 


Guidance 


23.2 


83.64 


62.68 


Warhead 


31.25 


28.26 


89 .91 


Sustainer 


87.51 


618.92 


149.29 


Booster 


26.31 


159.46 


206.20 


Motor Casing 


113.82 


59.50 


162.44 


Wings 




43.42 


211.40 


Canards 




9.88 


44.08 


Structure 


232.00 


172.38 


116.00 
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The complete missile weight is then 1264.62 lbs. 

The center of gravity from the above table is 

x _ = 139.44 in 
CG 

The center of gravity at the end of boost is 

> 

(x Q ) = 129.85 in 

G. DESIGN EXAMPLE (FINAL ANALYSIS) 

From the design parameters so far, a complete description 
of the missile can be determined. Figure (9-27) is a drawing 
of the missile. The launch conditions for the missile were 
specified as. 



“m = 1 * 5 

h = 10,000 ft 
W_ = 1264.62 

Cj 

x _ = 139.44 in 
cG 

From these conditions the input data for AEROl is shown in 
Table (9-III) . The format is the same as that of the output of 
AEROl and is printed as a check to ensure that the input data 
was entered properly. The output is shown in Table (9-IV) . 

Figure (9-28) is a plot of the coefficient of moment versus angle 
of attack for the launch condition and for the beginning of 
cruise. The center of gravity has moved forward approximately 
10 inches for the beginning of cruise so that the stability has 
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TABLE (9-IV) . AER01 Output 



\ 



O 



o 



n 

< 



n 

x 

> 





— 1 


O 


-H .O 




JO rH 


1*0 


CM 0 


n 


h* 


o 


ph ph 


rH 


J3 O 


O 


ph cr 


CJ 


O 


o 


0 

0 


O 


O O 


O 


O 0 



i 





<NI 


r-H 


(NO 


yo on 


f-H 


CM CM 


e_ 


jo 


O 


JO 


f-H 


yo 


-3* 3 


JO 


•JO 


'JO 


0 




















X 


0 


OO 


O *-H 


f-H 


f-H rH 


f-H 


i-H 


rH 




f-H 




r-H 


r-H 


f-H 


rH r— f 


f-H 


f-H 


f-H 


20 


0O 


O 


OO US 


30 


c- ph 


<2 


to 


Cl 


c_ 


ro 


O 


P-O 


O 


O 


ph cm 


O 


3 


CM 


0 


• 


















X 


yo 


O 

1 


(*0 


JO O 


O ph 


rH 


PH OO 


1- 


CM 


CM 


CM 


CM 


CM 


CM CM 


CM 


CM CM 


c. 


3 


3 


.3- 


3 


3 


3 3 


-3“ 


-3 


-3 


0 




















X 


-3- 


3 


3 


3 


3 


3 3 


-3- 


-3 


-3- 





to co 


CO co 


CO 


CO CO 


CO 


00 CO 






32 


Ph pH 


PH pH 


PH 


PH PH 


PH 


pH rn 






C. 


• • 
















O 


pH. pH 


PH PH 


pH 


rH ph 


PH 


pH rH 






X 


f-H f-H 


f-H r-H 


f-H 


rH rH 


rH 


rH f-H 








f-H rH 


r-H r-H 


JO 


3 cn 


3 


OO UO 








3 J T 


-3 -3 


-3 


’JO 0 


3 


CM rn 






< 


















CO 


0 

0 


0 

0 


O 


0 

0 


O 


1 — » rH 








UO O 


JO 0 


PH 


yo co 


f-H 


CM r-H 






31 


O — 


0 0 


3 


UO 0 


30 


ch 






0 


a 0 


0 *-< 


CM 


yo 3 


■JO 


PH 30 






c. 
















f-H 


t— 


0 0 


0 0 


O 


0 0 


O 


0 0 




3 


32 


0 0 


0 0 


O 


0 0 


O 


3 O 




O 


CO 
















3 




0 

0 


0 

0 


O 


0 0 


O 


0 

0 




• 


















O 


/ — 
















n 


r- 


0 0 


0 0 


O 


0 0 


O 


O O 




I— 


_J 


0 0 


0 0 


O 


0 0 


O 


O O 




rr 


CO 




















0 0 


0 0 


O 


0 0 


O 


3 O 




c 


















c 


a. 
















CJ 


>- 


rH 0 


pH UO 


-3 


CM CM 


rH 


O 3 








rH O 


r-H y> 


UO 


ph cn 


f-H 


P'0 3 






CO 


n 0 

t 


O O 


0 


a 0 


rH 


r-H rH 




CM 


















3 


G_ 


co 0 


CO O 


UO 


yo cm 


pro 


JO rn 




CM 


u 


r-H O 


r-H 


'.3 


CM 0 


O 


3 CC 




c 


















• 


CO 


0 

0 


0 

0 


O 


rH rH 


CM 


CM CM 




0 




1 














11 


















C_ 




CM O 


CM rH 


O 


CM rH 


O 


3 UO 




00 


— 


O O 


O rH 


VO 


O rH 


PH 


tC CO 




_ 


c 
















T~ 


CO 


0 

0 


O O 


0 


O rH 


rH 


cm yo 




c 


















CJ 


20 


3 0 


3 O 


PH 


O CM 


JO 


UO CM 








0 0 


O CM 


-3 


PH r-H 


UO 


0 3 






CO 














O 






0 

0 


O O 


0 

0 

1 


rH 


CM CM 


• 


c 




1 












O 


CM 
















II 


yo 


1— 


rH O 


PH JO 


CM O O 


yo 


3 O 


1- 


c 


23 


0 0 


O rH 


CM 


CM PPO 


3 


3 JO 


c 


• 


—1 














c: 


0 


CO 


0 0 


O O 


O 


O O 


O 


0 

0 


c_ 


N 




1 












c. 


h- 
















2J 


c 


* ■ 


PH O PH UO 


yo 


PM rH 


rH 


rH O 




c 


t- 


f-H O r-H -O 


UO Ph o 


rH p/0 JO 




CJ 


—1 


















a 


O O 
1 


O O 


0 


0 0 


r-H 


r-H r-H 








CM O 


-M -3 


PH 


CO O 


r-H 


CM PC 




0 


.23 


f-H O f-H CM 


yo 


3 O 




co 3 


O 


yo 
















• 


Ph 


CO 


0 

0 


O O 


O 


0 

0 


0 


0 

0 


O 


O 




1 












II 


• 
















1“ 


C 


CL 


CO O 


CO 3 


3 


CM Cl 


CM 3 O 


UL 


n 


72 


r-H O 


r-H 3 


JO 


CM O 


3 


3 CO 


r-' 


“2 


lo 














r- 


c 




O O 


O O 


3 


f-H f-H 


CM 


CM CM 


O 


c 




1 














CJ 


»- 


yo f-H 


yo cm 


f-H 


yo cm 


0 0 UO 






O 


JT -3 


3 UO *H 


0 JO 


CM 


rH CM 






1— 


















22 


0 

0 


0 

0 


O 


f-H rH 


CM 


yo 3 




UO 


O 
















PC 
















O 


0 


J- 


-3- O -3- Ph 


yo 


OO CM 


*0 


CM OO 


• 


p/0 


O 


0 0 0-3 


3 3 0 


CO 


r-H 3 


3 


• 


)- 














II 


C 


_J 


0 

0 


O rH 


"NJ 


yo 3 


5 

7 

8 


_J 


11 


CO 


1 














CO 
















— 


0 


_J 














c. 


c 


< 


CM O -Si -3 


0 


CO O 


Tl 


3 O 


CJ 


CJ 




1 






f-H 


rH 


f-H f-H 








<. f*-* 

h- v/0 




Z> cm 

C » 

>- h- 

CH 

era. 



264 



I 



o 



n 

< 



LU 

c: 



o 

o 



n 

2 

x 

> 





im. 


O 


fM m oi 


rH S -S' 


cD in 




rH 


O 


*H 


CM pH. 


3010 


rH OO 


2 














o 


O 


o 


O 


O O 


1 

2 

2 


m cm 




1 






1 1 


i i i 


1 1 




O CM 


o o o 


CO CM CM 


cn pa 


0_ 


>J\ 


-S 


•n 


o ^ 


CMf^rt 


CM I-H 


CJ 














X 


O 00 


o 


rH rH 


r-i r-1 


rH rH 




i-H 




H 


rH rH 


HH H 


rH rH 


DO 


CO O 


oo 


O r>- 


CM CO PA 


CD CD 


CL 


PA O 


pa 


rH rH 


CD -S' CD 


CM OO 


CJ 














X 


pn O 

1 


PA *A (D 


6 

7 

7 


OO OO 


t— 


CM 


CM 


CM 


CM CM 


CM CM CM 


CM CM 


Q_ 


-3* 


-S' 


-S' 


-S' -S' 


-3- -S' -S’ 


-S' -S 


CJ 














X 


-S' 


-S' 


-3- 


-S’ -S' 


-S' - s s 


-S -S' 





oo 


OO 


CO 


CO oo 


oo oo CO 


oo oo 




fM 


fM 


r-* 


p»- r>» 


r- r*» 


r»* r-*. 


51 














CJ 


p^. 


rM 


|M 


p- r>. 


psis Ps 


p^ 


X 


rH 


rH 


rH 


rH rH 


rH rH rH 


rH rH 




CD ID 


CD 


id 


PA CA OO 


cd in 




PA 


pa 




SIAN 


rH PA 


< 














<LJ 


O O 


o 


o o 


0 

0 

0 


rH CM 





tn o 


in 


r- o 


tn CM CM 


rH LA 






2 


n o 


tn 


CM rH 


OHfA 


P"*. rH 






<D 


















0 

0- 


o 


rH CM 


pa s tn 


CD CD 






Q 














CM 


P- 


o o 


o 


O O 


o o o 


O O 




-S' 


a 


o o 


o 


O O 


o o o 


O O 




CO 


o 














pa 




o o 


o 


O O 


0 

0 

0 


O O 




• 
















o 


c 














9 


1- 


o o 


o 


O O 


o o o 


O o 




H- 


-J 


o o 


o 


o o 


o o o 


O O 




CC 


CJ 


















0 

0 

0 


o o 


0 

0 

0 


0 

0 




C 
















CL 


2_ 














CJ 


p- 


-S' o 


-3* 


CD -S' 


O CD CM 


cd in 






z 


rH o 


rH 


CM S 


CD fM CD 


O CM 






CJ 


















O O 


o 


O O 


0 

0 

0 


rH rH 








1 












rH 
















O 


3- 


OO o 


OO 


CD 


rH CD <N 


O PM 




pa 


2 


rH O 


rH 


-S p>. 


rH S OO 


cm in 




c 


































CJ 


0 

0 


o 


0 

0 


1 

1 

1 


CM CM 




o 




1 












a 
















CL 




rH O 


rH 


oo pa 


CD rH tn 


OO CD 




cr. 


— 


O O 


O 


O CM 


s cd tn 


-s CO 




— 


o 
















o 


0 

0 


o 


0 

0 


0 

0 

1 


CM -S’ 




C 
















u 




CD O 


CD tn CD 


HtAPs 


|M OO 






-J 


O O 


o 


CM .S' 


oo cm rs 


-S’ rH 






o 












o 






O o 


o 


0 

0 


I 

I 

0 


CM -S’ 


• 


cn 




1 










o 


oo 














n 


CM 


p— 


-n o 


tn 


o tn 


o in cn 


-S’ oo 


I— 


o 


3 


o o 


o 


rH rH 


CM CM CM 


m pa 


c 


• 


— 1 












a 


c 


a 


o o 


o 


0 

0 


0 

0 

0 


o o 


Q_ 


9 




i 










Cl 


P- 
















o 


c_ 


-S' o 


S' 


CD -S’ 


o id m 


CD ID 




Cl 


p- 


»-» o 


rH 


CM -S 


CD CD 


o CM 




« 


















CJ 


o o 

» 


o 


0 

0 


0 

0 

0 


rH rH 








r- O 




tn cm 


CD rM m 


o r*- 




s 


D2 


o o 


o 


rH CM 


cm m s 


tn tn 


o 


CO 


-J 












• 


CO 


CJ 


o o 


o 


0 

0 


0 

0 

0 


0 

0 


o 


O 




1 










n 


• 














P- 


o 


h 


00 o 


OOtAN 


o tn cm 


O oo 


LL 


n 




rH O 


rH 


-3* a- 


H S OO 


cm tn 


< 


— fc 


Ij 












c: 


c 


L> 


o o 


O 


0 

0 


1 

1 

1 


CM CM 


<l 


Cl 




1 












CJ 


1- 


00 ID 


oo 


•n cd 


CD CO rH 


S CM 






O 


pa pa 


pa 


-S’ n 


00 CM CD 


OO O 






f— 
















3 


o o 


3 


0 

0 


I 

I 

0 


cm m 




oo 


o 














oo 














o 


CD 


1- 


pa o 


pa 


_S ID 


o in s 


O CD 


• 


CM 


o 


m o 


tn cm o 


O o CM 


CD CD 


o 


• 


t— 












II 


c 


-j 


o o 


o 


rH CM 


3 

k 

5 


CD OO 


-J 


■ 


CJ 


i 












cr 














— 


O 


-J 












a 


Cl 


< 


CM O 


CM 


-s id 


CO O CM 


-S’ CD 


CJ 


o 




» 






rH rH 


rH rH 







265 




266 






Figure (9-28). Static stability. 



increased as indicated in this figure. The performance objective 
from the beginning of the design was to provide a 31 g maneuver 
capability. Table (9- V) shows the output of AER01 for Mach 
numbers of 1.5 to 3.0 and control deflections of 0.0 to 30 degrees. 
From this output the trimmed normal force can be found and the 
corresponding load factor is then, 



'N, 



q s 



n = 



TR 



ref 



W 



The following values of maneuver load factor were found using 
AERO 1 . 
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deg. 



As indicated in the above table, the missile meets the maneuvering 
specifications of the operational requirements. 
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X. CONCLUSIONS AND RECOMMENDATIONS 



This thesis has presented the methods and general procedures 
for the conceptual design of tactical missiles. As mentioned, 
this is not necessarily the best design. Now that the design 
procedure has been attempted once, the reason for this is readily 
seen. The process is one of continuous compromise. As was shown 
in Chapter 7, the optimum wing for lift is not the best wing 
for minimizing drag. We also saw in Chapter 8 that increasing 
chamber pressure increases thrust at the expense of increased 
weight. This thesis has tried to point out some of these areas 
of compromise and present methods to deal with them. 

Areas which were not covered which need to be investigated 
in conceptual design are structures, radar cross-section and cost. 
With the increased emphasis on survivability and the decreasing 
budget, it becomes increasingly important to define the effects 
of these areas on design early in the process. 

The complexity of the design process and the need to obtain 
timely and accurate information have made it ideally suited for 
the digital computer. The AEROCF program used in this thesis 
is part of a large scale computer program (MISSYN) which con- 
sists of modules for each section of the design analysis. 

With a good understanding of the theory and methods used in 
missile design, the computer aided design program with graphics 
capability gives the designer the capacity to make intelligent 
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design interations and almost instantaneously see the effect 
of the change on all areas of design. A limited example of 
this can be seen in the use of AER01. A change in the per- 
formance requirements for maneuver capability would require 
a redesign of the lifting surfaces. A change in the lifting 
surface design would change the drag characteristics of the 
missile and therefore, the propulsion requirements. The AER01 
program coupled with a similar propulsion module would allow 
the designer to make the changes and instantly see the penalty 
or savings in propellant weight. 

One pass at the design has been accomplished in this thesis. 
As was seen throughout the process, decisions in one area 
affect the design in others. For this reason the design 
process becomes an iterative one. The final design of the 
first iteration is>.the baseline missile for the second itera- 
tion and the design is started again. By making several 
passes through the loop, the solution converges on the final 
design. 
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